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Abstract 

Free  space  optical  communication  (FSOC)  can  provide  the  modern  warfighter 
with  timely,  detailed  information  faster  and  more  securely  than  traditional  radio¬ 
frequency  (RF)  communications  systems.  However,  FSOC  systems  are  vulnerable  to 
unique  environmental  factors  and  design  constraints. 

The  propagation  of  an  FSOC  signal  through  atmospheric  turbulence  experiences 
random  fluctuations  in  intensity,  including  signal  fades  which  negatively  impact  the 
communications  link  performance.  This  research  develops  an  analytical  probability 
density  function  (PDF)  to  model  the  best  case  scenario  of  using  multiple  independent 
beams  to  reduce  the  intensity  fluctuations.  The  PDF  was  further  developed  to  account 
for  partially  correlated  beams,  such  as  would  be  experienced  by  beams  having  finite 
separation. 

The  PDF  was  validated  with  results  obtained  from  digital  simulations  as  well 
as  lab  experiments.  The  research  showed  that  as  the  number  of  transmitted  beams 
increases  the  probability  of  fade  decreases.  While  fade  probability  is  reduced  by 
adding  more  beams,  using  more  than  four  transmitters  does  little  to  improve  the 
overall  performance. 

Pulse  position  modulation  (PPM)  provided  significant  improvement  over  tradi¬ 
tional  fixed  threshold  on/off  keying  (OOK)  with  the  impact  of  signal  fading  reduced. 
When  integrating  the  use  of  PPM  with  four  transmitted  beams  in  a  simulated  airborne 
FSOC  system,  the  bit  error  rate  (BER)  was  reduced  by  a  factor  of  over  1  x  10-6  when 
compared  to  the  baseline  of  a  fixed  threshold  OOK  system  with  a  single  transmitted 
beam. 
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INTEGRATED  APPROACH  TO  FREE  SPACE  OPTICAL  COMMUNICATIONS 

IN  STRONG  TURBULENCE 

I.  Introduction 

Digital  high-speed  communication  is  a  vital  component  of  the  modern  battlefield, 
helping  to  provide  warfighters  with  timely,  detailed  information.  The  increase  in 
the  number  of  airborne  and  space-based  sensors,  as  well  as  their  increased  resolution 
and  hyperspectral  capabilities,  drives  the  need  for  an  expanded  wireless  communi¬ 
cation  capacity.  Optical  frequencies  have  the  potential  of  significantly  higher  data 
rates  than  traditional  RF  techniques,  but  optical  fibers  are  not  always  practical.  Free 
space  optical  communication  (FSOC)  systems  provide  the  opportunity  to  take  ad¬ 
vantage  of  the  higher  data  rates,  as  well  as  providing  a  lower  probability  of  intercept, 
while  avoiding  some  of  the  spectrum  regulation  issues  associated  with  radio- frequency 
(RF)  communication.  The  research  objective  is  to  further  develop  recent  methods 
for  mitigating  atmosphere-induced  signal  fluctuations  that  reduce  the  performance 
of  aerospace-based  FSOC.  Additionally,  the  research  further  develops  methods  to 
mitigate  the  effect  of  those  fluctuations  when  they  occur. 

When  laser  beams  are  transmitted  over  long,  turbulent  paths,  the  variation  in 
index  of  refraction  along  the  path  causes  the  beam  to  randomly  wander  and  scintillate, 
resulting  in  the  fluctuation  of  power  at  the  receiver  end  of  the  propagation  path.  At 
times  the  fluctuations  can  result  in  signal  fades,  causing  the  received  power  to  drop 
below  the  detection  threshold,  significantly  increasing  the  number  of  bit  errors.  For 
certain  engagement  geometries,  particularly  long  horizontal  paths,  traditional  phase- 
only  adaptive  optics  (AO)  cannot  adequately  correct  for  the  scintillation.  Other  non- 
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AO  approaches,  such  as  the  propagation  of  partially  coherent  beams  [5,65],  multiple 
receivers  [19],  and  polarization  diversity  [21]  have  also  been  suggested.  Additional 
methods  of  dealing  with  the  signal  loss  have  included  using  keying  schemes  based  on 
polarization,  pulse  positioning,  or  other  modulation  techniques  [35]. 

The  research  used  multiple  transmitters  separated  such  that  each  scintillates  in¬ 
dependently.  Louthain  and  Schmidt  showed  that  when  two  beams  are  separated  by 
a  distance  of  twice  the  Fresnel  zone  and  propagated  through  a  long  horizontal  path, 
they  scintillate  independently  [31,33].  When  the  two  beams  arrive  a  detector,  the 
effects  of  scintillation  average  out,  reducing  the  probability  of  signal  fades  and  their 
duration.  In  the  new  research  reported  here,  the  effects  of  additional  transmitters 
were  investigated  and  characterized.  For  practical  reasons,  the  transmitters  were  ar¬ 
ranged  to  minimize  the  amount  of  area  occupied.  Additionally,  to  mitigate  the  impact 
of  intensity  fluctuations,  various  modulation  techniques  were  investigated. 

To  demonstrate  the  improved  performance,  first  an  analytical  model  for  the  turbu¬ 
lence  irradiance  fluctuations  from  multiple  independent  beams  was  developed.  Build¬ 
ing  upon  that,  a  model  for  multiple  partially  correlated  beams  was  developed.  This 
model  was  then  combined  with  bit  error  rate  (BER)  models  to  predict  the  FSOC  sys¬ 
tem  performance.  These  models  were  then  compared  against  results  obtained  from 
digital  simulation.  Finally,  the  model  was  compared  with  results  obtained  from  a  lab 
experiment  simulating  a  FSOC  system  engagement. 

Key  contributions  of  this  research  include: 

•  An  analytical  model  for  the  probability  distribution  of  integrated  irradiance 
from  the  propagation  of  multiple  Gaussian  beams  through  a  turbulent  atmo¬ 
sphere 

•  Demonstration  of  improved  FSOC  system  performance  when  multiple  illumina¬ 
tors  are  combined  with  an  appropriate  modulation  scheme 
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•  Demonstration  of  decreased  integrated  irradiance  variance  for  multiple  illumi¬ 


nators  using  experimental  results,  verifying  the  analytical  model. 
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II.  Background 


This  research  builds  upon  previous  research  and  uses  principles  from  several  dis¬ 
ciplines,  including  digital  communications,  optical  propagation,  and  atmospheric  tur¬ 
bulence.  This  chapter  reviews  the  relevant  topics  in  each  area  with  application  to 
FSOC. 

2.1  Digital  communications 

A  model  of  a  basic  digital  communications  system  is  shown  in  Fig.  1.  The  com¬ 
munication  of  a  digital  bit  pattern  generally  begins  with  a  modulator.  The  modulator 
converts  some  number  of  bits  from  the  input  data  sequence  b%  into  a  physically  mea¬ 
surable  real-world  symbol  The  symbols  propagate  through  the  transmission 

channel  and  experience  the  channel  effects  which  typically  include  AWGN.  In  the 
case  of  FSOC,  channel  effects  also  includes  scattering,  absorption,  and  the  effects  of 
atmospheric  turbulence.  The  noise  added  within  the  receiver  is  also  included  in  the 
final  signal  r  (t).  From  the  received  signal,  a  demodulator  generates  a  bit  estimate 
that  is  passed  up  to  the  next  layer,  where  error  checking  can  occur. 

The  simplest  FSOC  modulation  scheme  is  the  simple  on/off  keying  (OOK)  tech¬ 
nique  where  the  laser  is  on  (signal  high)  to  transmit  a  symbol  corresponding  to  a  ‘1’ 
and  off  (signal  low)  to  transmit  a  symbol  corresponding  to  a  ‘O’.  There  are  then  four 
possible  scenarios  for  a  laser  receiver:  (1)  a  signal  high  was  sent  but  a  signal  low  was 


Figure  1.  Block  diagram  of  the  physical  layer. 
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estimated,  (2)  a  signal  high  was  sent  and  a  signal  high  was  estimated,  (3)  a  signal 
low  was  sent  and  a  signal  low  was  estimated,  and  (4)  a  signal  low  was  sent  but  a 
signal  high  was  estimated.  Scenarios  (1)  and  (4)  result  in  bit  errors.  The  BER  is 
typically  used  in  communications  systems  to  describe  system  performance,  with  the 
minimization  of  BER  being  the  goal.  The  optical  receiver’s  output  voltage  r  (t)  is 
modeled  by 

r(t)  —  I(t)h(t)G +  n(t),  (1) 

where  I  (t)  is  the  transmitted  irradiance  pattern  containing  the  communication  signal 
as  a  function  of  time,  h  (t)  is  the  channel  gain  as  a  function  of  time,  G  is  the  detector 
conversion  gain,  and  n  (t)  is  the  additive  noise  as  a  function  of  time.  The  channel  gain 
h  ( t )  is  a  random  process  due  to  the  atmospheric  turbulence.  The  additive  noise  n  (t) 
contains  the  aggregate  effects  of  shot  noise,  background,  and  receiver  electronics  which 
together  are  modelled  as  AWGN.  The  receiver  output  is  multiplied  by  a  reference 
signal  and  integrated  over  the  duration  of  symbol  to  generate  a  test  statistic  such 
that 

z=  f  sref  (t)  r  (t)  d  t,  (2) 

Jo 

where  z  is  the  resulting  test  statistic  in  units  of  Volts2  and  sref  (t)  is  the  reference 
signal.  For  OOK  the  reference  signal  is  simply  a  constant  value  over  the  symbol 
duration. 

A  maximum  likelihood  estimator  establishes  a  threshold  based  on  the  test  statistic 
probability  density  functions  for  the  transmission  of  a  ‘1’  p  (z\Hi)  and  the  transmission 
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of  a  ‘0’  p  (z\H0).  The  ratio  test  is  then 


A(*)  = 


p(z\Hi)  <  p(H0) 

p(z\H0)  >  p{Hi) 

Hn 


(3) 


which  is  used  to  determined  the  most  likely  symbol  sent.  Since  in  most  cases  the 
probability  of  a  ‘1’  or  ‘0’  is  equal  p  (H0)  =  p  (H i)  then  the  equation  simplifies  to 


if  p  (z\Hi)  >  p  (z\H0)  pick  Hi, 


(4) 


or 


if  p  (z\H0)  >  p  (z\Hi)  pick  H0. 


(5) 


For  binary  signals  the  test  statistic  is  compared  with  a  threshold  such  that  Eqs.  (4) 
and  (5)  simplify  to 

if  z  >  zt  pick  Hi,  (6) 


or 


if  z  <  zt  pick  Hq 


(7) 


The  optimal  threshold  for  OOK  is  determined  by 


*  =  ^ f  +  2  K  -  aS) !n  (  —  ) ,  (8) 


O’ l  Oq  <7 1  <7, 


1  u0 


O0 


where  po,i  is  the  mean  value  for  z  for  a  ‘O’  and  ‘1’  respectively  and  <r0)i  is  the  variance 
of  z  associated  with  a  ‘O’  and  T’  respectively  [6, 14,30]. 
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2.2  Free  space  optical  propagation 


For  FSOC,  light  is  used  to  create  the  symbol,  representing  a  bit  or  packet  of  bits, 
which  is  propagated  through  an  open  medium.  Before  considering  the  turbulent  ef¬ 
fects  that  the  medium  has  on  the  propagated  light,  vacuum  propagation  is  considered. 
Since  most  laser  beams  can  be  well-represented  by  Gaussian  beams,  the  propagation 
of  Gaussian  beams  is  also  discussed. 


2.2.1  Vacuum  propagation. 

If  a  propagation  medium  is  linear,  isotropic,  homogeneous,  non-dispersive,  and 
nonmagnetic,  then  Maxwell’s  equations  can  be  reduced  to  the  scalar  wave  equation 


V2w  (P,  t) 


n2  d2U  (P,  t) 
c2  dt2 


0, 


(9) 


where  n  is  the  index  of  refraction  of  the  medium,  c  is  the  speed  of  light,  u  represents 
any  of  the  scalar  field  components  at  a  position  P  and  a  time  t  [20].  If  it  is  further 
assumed  that  the  propagated  light  is  monochromatic  then  the  scalar  field  can  be 
described  as 

u  (P,  t)  =  A  (P)  cos  [27 TVt  +  <f>  (P)] ,  (10) 

where  A  (P)  and  (f)  (P)  are  the  amplitude  and  phase  of  the  wave  at  a  position  P,  and 
v  is  the  optical  frequency.  More  compactly,  Eq.  (10)  becomes 


u  (P,  t)  =  Re  {U  (P)  exp  (— f27rz/t)}  , 


(11) 


where  Re{}  signifies  the  “real  part”,  and  U  (P)  is  the  complex  pliasor 


U(P)  =  A  (P)  exp  [-i(j>  (P)] . 


(12) 
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Since  the  time  dependence  is  always  exp  (—i2iri't)  for  a  traveling  wave,  then  U  must 
obey  the  time-independent  Helmholtz  equation 


(V2  +  k2)  U  =  0, 


(13) 


where  k  =  2tt/\  with  A  being  the  field  wavelength.  Solving  Eq.  (13)  using  Green’s 
theorem  and  the  free-space  Green’s  function  provides  the  Huygens-Fresnel  princi¬ 
ple  [20] 

U(Pl)  =  \  !>«*,«.,  (14) 

*  A  JJ  roi 

s 

where  U  (Pq)  and  U  (Pi)  are  the  incident  and  diffracted  fields,  6  is  the  angle  between 
the  aperture  normal,  the  vector  r*oi  pointing  from  Pq  to  Pi,  and  £  defines  the  aperture 
area. 

The  Huygens-Fresnel  principle  can  be  re-written  in  Cartesian  coordinates  as 


u^=kSfu 


(15) 


by  noting  that  the  term  cos  6  is  given  by 


cos  6 


z 

r0 1 


(16) 


and  the  distance  r0i  is  given  by 

roi  =  \Jz2  +  (x  -  O2  +  (y  -  v)2-  (17) 

Approximations  include  scalar  diffraction  and  that  r0i  >  A.  Further,  by  taking  the 
first  two  terms  of  the  binomial  expansion  of  Eq.  (17),  the  resulting  Fresnel  diffraction 


integral  can  be  expressed  as 


U{x,y) 


g ikz 

iXz 


U  R,  >)) 


(18) 


where  U  (£,  rj)  is  the  original  held  and  U  (x,  y )  is  the  result  of  propagating  a  distance 
of  z.  Re-writing  in  terms  of  a  Fourier  transform  and  performing  some  algebra  yields 


Akz 


u 


(x'v)=j\/ 


( x2+y 2) 


FW&rDe'*- 


{ew 


(19) 


This  form  shows  that  the  propagated  held  is  merely  a  scaled  Fourier  of  the  original 
held  multiplied  by  a  quadratic  phase  factor  demonstrating  that  the  Fresnel  diffraction 
integral  is  useful  for  linear  systems  analysis. 


2.2.2  Propagation  of  Gaussian  beams. 

Most  optical  beams  propagating  in  free  space  can  be  modeled  as  purely  transverse 
electro-magnetic  (TEM)  modes  [63].  The  lowest-order  Gaussian  beam  mode,  TEM0o, 
is  typically  used  to  model  a  laser  beam,  with  limiting  cases  including  plane-waves  and 
spherical-waves.  The  Gaussian  beam  mode  is  found  by  solving  the  wave  equation. 
This  section  describes  the  vacuum  propagation  of  a  TEMoo  mode  Gaussian  beam. 

A  Gaussian  laser  can  be  parameterized  by  its  1/e  radius  Wo,  its  radius  of  curvature 
F0,  and  its  wavelength  A.  The  optical  held  at  the  source  plane  of  the  TEM00  wave 
can  then  be  dehned  as  [1] 

/  r2  ikr2\  (  1  \ 

U0  ( r,z  =  0)  =  a0  exp  f -  Tjjr  J  =  aoexp  \  --a0kr2j  ,  (20) 

where  r  =  (x2  +  y 2)1^  is  the  radial  distance  from  the  beam  center,  z  is  the  propaga¬ 
tion  distance,  ao  is  the  optical  held  amplitude  in  (W/m2)  ,  and  «o  is  the  complex 
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beam  parameter  related  to  spot  size  and  phase  front  radius  such  that 


(21) 


The  propagation  beam  parameter  p  (z)  is  expressed  in  the  form 


p  {z)  =  1  +  ia0z  =  @o  +  1A0, 


(22) 


where  ©o  and  Ao  are  the  input  beam  parameters  defined  as  the  real  and  imaginary 
parts  of  p  (z)  such  that 


The  parameter  @0  is  also  known  as  the  curvature  parameter  while  A0  is  the  Fresnel 
ratio  at  the  input  plane.  For  a  fixed  path  length  and  initial  radius  of  curvature  Fq, 
the  curvature  parameter  for  collimated,  convergent,  and  divergent  beams  are  @0  =  1, 
@o  <  1,  and  @o  >  1,  respectively.  The  limiting  cases  of  a  plane- wave  and  spherical- 
wave  (i.e.  point  source)  are  modeled  with  Wo  =  exam  and  Wq  =  Orn  resulting  in 
A0  =  0  and  A0  =  oo,  respectively. 

The  propagated  Gaussian  beam,  determined  by  solving  the  wave  equation,  is 


(24) 


where  p,  W,  and  F  represent  the  longitudinal  phase  shift,  spot  size  radius,  and  radius 
of  curvature  at  a  position  0  along  the  propagation  path.  These  are  defined  in  terms 
of  the  input  beam  parameters  as 


ip  =  tan  1 


(25) 
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(26) 


w  =  was/e$ -t- a$, 

Fa  (9?  +  Aq)  (0o  -  1) 


®0  +  Aq  —  00 

From  this,  the  output-plane  beam  parameters  are  defined  as 


(27) 


^  ,  2  2  z 

+  F’  ~  kw 2' 


(28) 


2.3  Atmospheric  turbulence 

Atmospheric  turbulence  degrades  the  FSOC  performance  in  several  ways.  The 
propagated  beam  is  dispersed  and  moved  about,  limiting  the  amount  of  light  incident 
on  the  receiver.  Additionally,  higher-order  intensity  fluctuations  vary  the  amount  of 
light  incident  on  the  receiver.  Ultimately,  the  turbulence  causes  fades  in  the  signal 
intensity  resulting  in  errors  in  the  bit  estimates.  Even  seemingly  short  fades,  like 
lmsec,  for  a  lGbit/sec  system  could  result  in  as  many  as  500,000  bit  errors. 


2.3.1  Atmospheric  fluctuations  of  index  of  refraction. 

Optical  turbulence  in  the  atmosphere  results  from  fluctuations  in  the  index  of 
refraction  due  to  turbulent  air  motion  through  the  variation  of  the  pressure,  temper¬ 
ature,  and  velocity  of  air.  Light  propagating  through  this  turbulent  medium  expe¬ 
riences  distortion,  both  spatially  and  temporally,  resulting  in  the  twinkling  of  stars, 
loss  of  resolution  for  imaging  systems,  and  the  distortion  of  propagating  laser  beams. 

The  ultimate  source  of  the  turbulence  arises  from  the  heating  and  cooling  of  the 
surface  of  the  earth  caused  by  sunlight  and  the  diurnal  cycle,  which  cause  large-scale 
variations  in  air  temperature.  Kolmogorov  developed  a  model  for  the  velocity  struc¬ 
ture  function  and  power  spectral  density  (PSD)  of  turbulent  velocity.  He  suggested 
that  the  kinetic  energy  in  large-scale  motion  gets  transferred  to  smaller  scale  motions 
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resulting  in  a  transition  from  turbulent  to  laminar  flow  [47].  This  allowed  others  to 
apply  similar  analysis  to  temperature  and  thus  derive  a  model  of  the  refractive-index 
structure  function  and  PSD  [12,37]. 

At  optical  wavelengths,  the  dependence  of  the  index  of  refraction  of  air  on  tem¬ 
perature  and  pressure  is 

77  6  P 

n  =  n0-\ - —  x  1(T6,  (29) 

where  no  =  1,  T  is  the  air  temperature  in  Kelvins,  and  P  is  the  air  pressure  in 
millibars.  The  temperature  fluctuations  normally  dominate  the  index  of  refraction 
variations.  However,  in  boundary  layer  regions,  such  as  around  an  optical  window  on 
an  airframe,  the  pressure  variations  can  be  large  enough  to  cause  significant  changes 
in  the  index  of  refraction. 

The  statistical  distribution  of  the  size  and  number  of  the  turbulence  eddies,  regions 
of  uniform  index  of  refraction,  is  characterized  by  the  power  spectral  density  of  n  (r), 
denoted  by 

3  /»oo 

/  Bn  (r)  exp  (— in  ■  r )  dr,  (30) 

J  —  OO 

where  Bn  (r)  is  the  autocorrelation  of  n  (r)  and  r  is  a  position  vector.  The  indepen¬ 
dent  variable  k,  is  the  spatial  wavenumber  vector  with  orthogonal  components  along 
the  x,  y,  and  2  directions.  Under  conditions  where  homogeneous,  isotropic  turbu¬ 
lence  is  assumed,  the  power  spectral  density  of  the  index  of  refraction  fluctuations  is 
a  function  of  the  scalar  wavenumber 

«=  \J  +  k2z, 

which  is  related  to  the  isotropic  scale  size  by  /  =  2i t/k. 

Standard  Kolmogorov  theory  does  not  predict  a  mathematical  form  for  (n) 
outside  of  the  region  where  larger  eddies  break  up  into  smaller  ones.  The  quantity  Lq 


(31) 
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is  referred  to  as  the  outer  scale  and  /0  the  inner  scale,  between  which  this  assumption 
holds  true.  When  k  is  within  the  range  2tt/L0  -C  k  <C  27 r//0,  the  form  of  $„(/«) 
predicted  by  Kolmogorov  theory  is  [47] 


(k)  =  0.0 33C2„K-u/3, 


(32) 


where  U2  is  the  index  of  refraction  structure  constant  with  units  of  m  2//3  and  the 


superscript  K  is  used  to  denote  the  Kolmogorov  spectrum.  When  inner  and  outer 
scales  are  considered,  when  operating  outside  the  range  where  2n/L0  <C  k  <C  27 r/Z0, 
the  modified  von  Karman  spectrum  model  is  used  and  is  represented  by  [47] 


(33) 


where  kq  =  2tt/L0,  nm  =  5.92 /l0,  and  the  superscript  V  denotes  the  von  Karman 
spectrum. 

For  most  problems  of  interest,  the  atmospheric  turbulence  strength  varies  along 
the  propagation  path.  In  this  case,  the  index  of  refraction  structure  parameter  be¬ 
comes  C\  (z)  where  z  is  longitudinal  position  along  the  path.  Using  this  notation, 
the  Kolmogorov  and  von  Karman  spectrums  from  Eqs.  (32)  and  (33)  become 


**  (ic,z)  =  0.033Cj(s)  K  11/;:, 


(34) 


and 


(35) 


respectively. 

A  commonly  used  model  of  (72  is  the  Hufnagcl- Valley  turbulence  prohle  described 
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by  [47] 

C2n(h)  =  5.94  x  10-53  (w/27)2  h10  exp 

+  2.7xlO-1«exp(-IA5)+^exp(-A),  (36) 

where  h  is  the  height  above  the  ground  in  meters,  A  is  a  free  parameter  that  sets  the 
turbulence  strength  near  the  ground  with  units  of  m-2/3,  and  v  is  a  free  parameter 
that  sets  the  high-altitude  wind  speed  with  units  of  m/s.  Commonly  used  values  are 
A  =  1.7  x  10-14m-2/3  and  v  =  21m/s,  often  referred  to  as  the  HV5/7  model.  Other 
models  include  the  SLC-Day  model,  the  Modified  Hufnagel-Valley  model,  and  the 
Greenwood  model.  A  comparison  between  the  models  is  shown  in  Fig.  2. 
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2.3.2  Atmospheric  turbulence  parameters. 


The  effects  of  atmospheric  turbulence  are  observed  through  several  optical  phe¬ 
nomena.  These  include  phase  and  amplitude  effects,  which  cause  laser  beams  to 
spread,  wander,  and  scintillate.  The  quantities  used  to  describe  the  accumulated 
phase  effects  include  the  atmospheric  coherence  length  and  the  isoplanatic  angle, 
while  quantities  describing  the  amplitude  effects  include  log-amplitude  variance  and 
intensity  variance. 

The  development  of  the  phase  and  amplitude  effects  starts  with  the  stochastic 
Helmholtz  equation 

V2f/  +  k2n2  (r)  U  =  0,  (37) 

where  the  refractive  index  is  expressed  as  a  function  of  position  r 

n  ( r )  =  n0  +  ni  (r) ,  (38) 

and  Hq  =  ( n(r ))  =  1.  Using  the  Rytov  method,  one  of  several  methods  to  solve  for 
the  turbulence  effects,  the  field  of  the  electromagnetic  wave  is  written  as 

U  (R)  =  U  (r,  L)  =  U0  (r,  L)  exp  (r,  L)] .  (39) 

The  complex  phase  perturbation  due  to  the  turbulence,  '0?  is  of  the  form 


^  (r,  L)  =  0i  (r,  L)  +  02  (r,  L)  +  •  ■  ■  ,  (40) 

where  0i  (r,  L)  and  02  (r,  L)  represent  the  first-order  and  second-order  perturbations, 
respectively.  Various  statistical  moments  of  the  perturbations  can  then  be  computed, 
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for  example 


£1  (0,0) 


(^2  (r,L))  +  -  (^2(r,L)) 


r»L  /»oo 


—  27 T  k 


2  1,2 


(k,  z)  dftdz, 


'0  Jo 


(41) 

(42) 


where  then  the  appropriate  spectrum  (Kolmogorov,  von  Karman,  etc.)  could  be  used. 
This  moment  and  others  then  can  provide  the  basis  for  parameters  used  to  describe 
the  phase  and  amplitude  effects  of  propagation  through  turbulence. 

The  effect  of  beam  spreading  is  related  to  the  limiting  effect  of  the  atmosphere 
on  imaging  resolution.  Without  atmospheric  turbulence,  the  resolution  of  an  imaging 
system  is  determined  by  the  aperture  diameter,  known  as  the  diffraction  limit.  When 
atmospheric  turbulence  is  considered,  imaging  system  resolution  can  be  approximated 
by  replacing  the  aperture  diameter  with  the  atmospheric  coherence  diameter  ro  also 
known  as  the  Fried  parameter  [16].  Using  the  Kolmogorov  spectrum,  for  a  plane- wave, 
this  results  in  [49] 

r  1  3/5 


T*  0  ,pw 


0.423fc2  /  C2(z)ck 


(43) 


and  for  a  spherical- wave  the  result  is 


^0,sw 


(44) 


where  k  =  2ti/X  is  the  wavenumber  and  L  is  the  propagation  path  length.  The 
plane-wave  and  spherical-wave  results  provide  limiting  cases  for  a  Gaussian  beam. 
Therefore,  it  is  common  to  use  the  Gaussian  beam  parameters  to  determine  whether  it 
behaves  more  like  a  planar  or  spherical-wave  and  then  use  the  corresponding  equation. 

Compensation  for  the  phase  distortions  in  the  atmosphere  relies  upon  an  idealized 
point  source  beacon  in  the  vicinity  of  the  object  to  be  observed.  How  close  this  point 
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source  needs  to  be  is  determined  by  the  turbulence  along  the  path.  The  separation 
is  given  as  the  isoplanatic  angle  and  is  found  via  [49] 


e  0  = 


2.91  k2 


(45) 


As  long  as  the  reference  and  the  object  to  be  observed  are  within  this  angular  ex¬ 
tent,  phase  compensation  based  on  the  reference  can  effectively  compensate  for  the 
distortions  that  the  image  experiences. 

Atmospheric  distortion  of  a  laser  beam  also  includes  amplitude  effects,  also  known 
as  scintillation.  As  the  beam  propagates  through  the  atmosphere,  the  cumulative  ef¬ 
fects  of  phase  aberrations  result  in  constructive  and  destructive  interference  at  the  re¬ 
ceiver,  which  are  observed  as  intensity  variations.  For  laser  communications  systems, 
these  effects  are  particularly  important  since  they  can  cause  signal  fades  resulting 
in  bit  errors.  This  phenomenon  can  be  characterized  by  the  log-amplitude  variance 
cr£,  commonly  referred  to  as  the  Rytov  number  1Z  and  is  calculated  using  the  Rytov 
approximation.  Turbulence  strength  is  often  expressed  in  terms  of  the  Rytov  number. 
The  plane-wave  log-amplitude  variance  at  the  receiver  is  [1] 

=  0.5631  k7^  [L  C2n  (z)  (L  -  zf 6  dz,  (46) 

•vo 


and  the  spherical-wave  log-amplitude  variance  is 

v2x,Pw  =  0.5631  k7/6  Cl  (z)  (-0  7  (L  -  z)5/ 6  dz.  (47) 

The  expressions  in  Eqs.  (46)  and  (47)  are  only  valid  for  weak  turbulence,  where 
the  log-amplitude  variances  is  <  0.25.  For  moderate  to  strong  turbulence,  when 
<7^  >  0.25,  the  scintillation  begins  to  saturate  and  the  Rytov  approximation  no  longer 
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holds.  For  horizontal  propagation  (i.e,  constant  turbulence)  of  plane-waves,  the  log- 
amplitude  variance  reduces  to 


cr 


2 

x,pw 


0.3071  k7'6C2nLn/6. 


(48) 


The  detector  at  the  receiver  measures  irradiance,  not  the  field  amplitude  directly. 
For  weak  turbulence,  the  normalized  irradiance  variance  (also  known  as  the  scintilla¬ 
tion  index)  erf  is  [1] 


(/2(r))  (/(r))2  (J2(r)) 

</(r))2  </(r)>2  (J(r))2 

=  exp  [erf  (r)]  -  1 


(49) 

(50) 


=  4 erf  (r)  for  erf  <  0.25 


(51) 


=  ofn /(r)  for  erf  <  0.25, 


(52) 


where  r  is  the  radial  distance  from  the  optical  axis  and  ofn/  (r)  is  the  log-irradiance 
variance.  Based  on  this  relationship,  the  plane-wave  scintillation  index  is 

ofn/  =  2.251  k7/6  f  (F2  (z)  ( L  -  z)b/&  &z,  (53) 

Jo 

and  for  a  spherical-wave  is 

”1,  =  2.25U:7/6  £  Cl  (z)  Q5'6  (I  -  zf6  (54) 

Andrews  and  Phillips  define  the  plane- wave  scintillation  index  for  a  constant  C2  profile 
as  the  Rytov  variance  erfj,  which  is  not  to  be  confused  with  the  Rytov  number  [1], 
From  Eq.  (53),  the  Rytov  variance  is  defined  as 


4  =  4  =  1. 23k7/6  C„Ln/6. 


(55) 
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For  strong  turbulence,  Rytov  numbers  larger  than  0.25,  the  phase  variance  continues 
to  increase  while  the  amplitude  variance  saturates.  Several  theories  have  been  used 
to  predict  this  behavior  as  the  Rytov  number  increases  [1].  Extended  Rytov  theory 
breaks  the  normalized  irradiance  into  two  factors  according  to 


(56) 


where  /  is  the  normalized  irradiance,  X  arises  from  large-scale  turbulent  eddy  ef¬ 


fects,  and  Y  arises  from  statistically  independent  small-scale  turbulent  eddy  effects. 
The  small-scale  and  large-scale  variances  are  typically  modeled  as  Gamma  random 
variables,  resulting  in  /  being  a  Gamma-Gamma  random  variable.  The  normalized 
((/)  =  1)  Gamma-Gamma  distribution’s  probability  distribution  function  (PDF)  is 
given  by 


(57) 


I  >  0, 


where  a  =  l/cr^,  (3  =  1  /crjj,  F  (■)  is  the  Gamma  function  and  Kp  (■)  is  the  modified 
Bessel  function  of  the  second  kind.  More  often,  the  cumulative  distribution  function 
(CDF)  is  used  to  determine  fade  statistics,  since  it  describes  the  probability  of  re¬ 
ceiving  less  than  a  given  threshold.  The  CDF  for  the  Gamma-Gamma  function  can 
be  determined  analytically  to  be 


o 

71 


sin  [7r  ( a  —  /3)\  T  (a)  T  ((3) 


(58) 
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where  It  is  a  threshold,  /3\  =  /3  —  a  +  1,  aq  =  a  —  (5+1  and  i F2  is  a  generalized 
hypergeometric  function. 

The  total  scintillation  index  then  takes  the  form  of 


a)  =  exp  (a -  1  =  exp  (a^x  +  irj;,,)  -  1, 


(59) 


where  crfnX  and  aj^y  are  the  variances  of  the  small-  and  large-scale  irradiance  fluc¬ 
tuations,  respectively.  As  with  the  equations  for  weak  turbulence,  the  equations  for 


01 


In  A' 


and  ofny  take  different  forms  depending  on  the  type  of  source.  For  Kolmogorov 


turbulence  (i.e,  zero  inner-scale  and  infinite  outer-scale)  the  log  variances  are  given 
by  [1] 


‘TnA  — 


0.49a! 


i  +  i.h42/5)7/6 


aR  ^ 


and 


a\nY  — 


0-51  a% 


1  +  0.69ai2/5) 


7/6 


for  a  plane-wave  and 


2 

°inX  — 


0.20a! 


l  +  0.19a!2/5)7/6 


and 


ainY  — 


0.20a! 


1  +  0.23a!2/5) 


7/6 


0.51a!, 

In  2, 

4» 

0.20a!, 

1.37 

4/5  > 

4» 

aR 

0.20  a!, 

4« 

In  2, 

<4» 

(60) 


(61) 


(62) 


(63) 


for  a  spherical- wave. 

For  practical  applications,  such  as  for  FSOC,  the  irradiance  is  averaged  across  the 
receiver  aperture,  rather  than  determined  at  a  single  point,  and  the  light  is  focused 
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onto  a  receiver  that  detects  the  total  power  in  the  aperture.  This  allows  the  receiver 
to  capture  more  energy  from  the  propagating  beam,  and  also  averages  out  the  spa¬ 
tially  varying  scintillation  effects.  For  a  Gaussian  beam,  and  assuming  Kolmogorov 
turbulence,  the  aperture-averaged  log-variance  of  x  is  [2] 


<A°)  -  0.49a?  (i-^  +  E6? 


X 


1  7/6 


1  +  0A0r]x  (2  —  ©i)  /  (Ai  +  Gg) 


(64) 


where 


Vx  = 


(|-l0i  +  le?)  6/Wi)12/7 


1  +  0.56^2/5) 


(65) 


The  aperture-averaged  log-variance  of  y  is 


a 


in  y 


(D)  = 


1-27  afriy' 


1  +  0A0r]y/  (Ai  +  Gg) 


Vy  < 


(66) 


where 

/  \  12/5 

Vy  =  3  (^rj  (l  +  0.69c42/5)  .  (67) 

The  variance  is  determined  via  a2  y  =  exp  —  1.  The  parameter  a2B  is  the  Rytov 

variance  for  a  beam  wave,  which  is  approximated  as 


a 


B 


^  3.86of{0.40  [(1  +  2©i)2  +  4A^] 


5/12 


X  COS 


5 

-tan 
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(68) 

(69) 


where  a2  is  the  plane-wave  Rytov  variance  given  as  erf  =  1.23C^fc7//6Ln/6.  The  values 
@i  and  Ai  are  the  curvature  parameter  and  Fresnel  ratio  at  the  receive  plane  for 
vacuum  propagation,  which  given  in  terms  of  their  respective  values  at  the  source 
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plane,  are  [1] 


(70) 


9i 


0q 

eg  +  A§ 


and 


Ai  = 


An 


eg  +  Ag' 


(71) 


The  parameter  @0  is  simply  defined  as  @0  =  1— ©o-  Finally,  the  value  f2G  characterizes 
the  finite  size  of  a  Gaussian  lens  and  is  given  to  be  Gg  =  2 L/IcWq  where  Wq  = 
(D2/8)1//"  with  D  representing  the  receive  aperture  diameter. 

The  temporal  nature  of  atmospheric  turbulence  is  typically  described  using  the 
Tyler  and  Greenwood  frequencies  given  by  [47] 


and 
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(72) 

(73) 


(74) 

(75) 


where  |V|  is  the  mean  wind  velocity  along  the  path.  The  majority  of  tilt  jitter 
fluctuates  at  rates  below  the  Tyler  frequency.  Similarly,  the  majority  of  higher- 
order  and  lower-order  phase  disturbances  fluctuates  at  rates  below  the  Greenwood 
frequency.  This  provides  a  design  point  for  tilt  and  higher-order  compensator’s  closed- 
loop  bandwidths.  The  wind  velocity  in  Eqs.  (72)  and  (74)  is  frequently  described  by 
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the  Button  wind  model  given  as  [1] 


V  ( h )  =  ujgh  +  14  +  30  exp 


(76) 


where  Vg  is  the  ground  wind  speed  and  ug  is  the  slew  rate  of  a  ground-based  telescope 


tracking  an  aerial  target. 

2.3.3  Turbulence  effects  on  FSOC. 

Beam  spreading  and  beam  wandering  alone,  without  considering  the  impact  of 
scintillation,  can  cause  the  beam  to  attenuate  and  wander  off  the  receiver  as  a  re¬ 
sult  of  atmospheric  turbulence.  When  this  occurs,  the  receiver  experiences  a  fade. 
For  horizontal  propagation  (constant  C%),  the  Gaussian-beam  wander  displacement 
assuming  infinite  outer  scale  is  [1] 


where  2^1  is  a  hypergeometric  function  and  ©o  is  the  curvature  parameter  defined  in 
Eq.  (23).  When  the  beam  is  collimated,  @0  =  1,  then  Eq.  (77)  reduces  to 


(rc)  =  2.42C^L3W0-,/3 


(78) 


and  when  the  beam  is  focused,  @0  =  0,  then  it  reduces  to 


(rc)  =  2.72  C3L3W0-1/3 


(79) 


If  the  turbulence  is  assumed  to  be  isotropic,  then  the  Cartesian  components  of  the 
displacement  can  be  treated  as  identically  distributed,  zero-mean,  Gaussian  random 
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variables  xc  and  yc.  They  are  related  to  the  displacement  via 


rc  =  \/x2c  +  y2c. 

(80) 

The  variances  of  xc  and  yc  can  then  be  found  by 

<c2>  =  {xl  +  vl), 

(81) 

=  <*2>  +  <2/2>, 

(82) 

=  <*2>  +  <*2>, 

(83) 

=  2<x2>, 

(84) 

since  (x2)  =  (?/2).  The  beam  displacement  follows  a  Rayleigh  distribution  with  a  PDF 
of 

^ 

/«  M  =  ^exp^— 
with  a  cumulative  distribution  function  (CDF)  of 

Fr  (r) = L  /b  K)  d? = ^ exp  (^) = 1  ~ exp  (sO  •  (86) 

where  a2  =  (x2)  =  (r/2)  =  1/2  (r2).  By  way  of  example,  a  collimated  Gaussian  beam 
propagating  100km  through  horizontal  turbulence  at  C*2  =  1  x  10-17m-2/3  with  an 
initial  beam  waist  of  Wo  =  5cm,  results  in  a  PDF  and  CDF  of  displacement  shown 

in  Fig.  3. 

When  combined  with  the  turbulence-induced  beam  wander,  short-term  beam  ra¬ 
dius  Wst  at  the  receiver  can  be  used  to  determine  the  average  irradiance  at  the 
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Figure  3.  PDF  and  CDF  for  5cm  beam  propagating  100km. 


receiver.  The  short-term  beam  radius  is  given  by  [1] 


WST  =  W 


\ 


1  +  1.33a2  A5/6 


1  -  0.66 


An 


1  +  An 


1/6' 


(87) 


for  a  collimated  beam  and  by 


WST  =  wJ\  +  0.35<r|A5/«. 


(88) 


for  a  focused  beam.  The  long  term  beam  radius  can  then  be  found  via  [1] 

WLT  =  sjwfr  +{r|),  (89) 

where  r2  is  found  from  Eq.  (78)  or  Eq.  (79).  Continuing  from  the  previous  example, 
the  short  term  beam  radius  is  Wst  —  1.005m,  the  mean-square  displacement  of  the 
beam’s  centroid  is  (r2)  =  0.0657m2,  and  the  long  term  beam  radius  is  Wlt  =  1.037m. 

The  received  power  of  each  possible  beam  location  is  found  by  integrating  the 
displaced  received  irradiance  over  a  circular  telescope  pupil.  The  resulting  captured 
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power  P  ( xc ,  yc)  for  a  beam  centered  on  the  Cartesian  coordinates  (xc,  yc)  is  given  by 


A 


(90) 


where  Jr  is  the  irradiance  at  the  receive  plane  and  A  represents  the  receive  aperture 
area.  Relating  the  received  power  at  a  given  position  and  the  probability  of  each 
position,  a  CDF  can  be  created.  Using  the  same  scenario  and  turbulence  conditions 
with  1W  of  transmitted  power  and  a  10cm  diameter  receiver,  the  resulting  CDF  is 
shown  in  Fig.  4. 

The  fade  probability  is  found  by  determining  the  probability  at  the  detection 
threshold.  This  calculation  ignores  the  higher-order  effect  of  scintillation,  which 
causes  spatial  variation  in  intensity  across  the  receive  plane.  Using  wave-optics  sim¬ 
ulation,  a  scintillation  included  CDF  can  be  calculated  via  numerous  independent 
random  realizations  of  the  atmosphere.  Using  the  same  scenario  and  turbulence  con¬ 
ditions  from  before,  the  resulting  CDF  is  shown  in  Fig.  5. 

In  addition  to  the  above  method  of  determining  the  long-term  turbulence-induced 
beam  waist,  Andrews  and  Phillips  give  the  effective  beam  width  We  as  a  result  of 
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Figure  4.  CDF  for  tilt-only  disturbance. 


Figure  5.  CDF  for  scintillation  included  disturbance. 
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propagating  through  turbulence  as  [2] 


We  =  w(l  +  1.63aJ2/5A^  V2  .  (91) 

If  all  of  the  beam  displacement  is  due  to  atmospheric  turbulence,  and  not  pointing 
errors  induced  at  the  transmitter,  then  the  effective  beam  waist  is  the  same  as  the 
long-term  beam  waist  Wlt  =  We.  This  formulation  proves  useful  to  calculate  the 
mean  irradiance  at  the  receiver  plane. 


2.3.4  Anisoplanatic  separation. 

In  previous  work,  Louthain  determined  the  analytic  log-amplitude  and  phase 
structure  functions  for  a  horizontal  path  starting  with  the  von  Karman  PSD 


(k,z) 


0.033C*  (z) 
{k2  +  Ko)11/6 


(92) 


where  k  is  the  three-dimensional  radial  spatial  frequency  and  k0  =  2n/L0  [1,47].  In 
general,  a  structure  function  associated  with  a  random  process  x  (t)  =  /ix  (t)  +  X\  (£) 
is  defined  by 

Dx  (ti,t2)  =  ([x  (ti)  -  x  (t2)]2)  ,  (93) 


where  fj,x  is  the  mean  value  of  x.  If  fj,x  is  slowly  varying,  then  the  difference  in  the 
means  can  be  ignored,  making  it  useful  for  the  characterization  of  random  processes 
with  stationary  increments  such  as  this.  When  the  process  is  stationary  in  the  wide 
sense,  then  the  structure  function  reduces  to 


Dx  ( t )  =  2  [i4  (0)  -  Rx  (t)]  , 


(94) 


where  Rx{r )  is  the  autocorrelation  of  x(t).  The  spatial  log-amplitude  structure 
function  is 


Dx(d) 


27tL°^(XLk2\]  KdK  (q^ 

\L^Sm{2nLl)\  («2  + ^ 


where  d  is  separation  the  separation  distance  between  points  of  interest  in  the  receive 
plane  and  Jq  is  the  zeroth-order  Bessel  function  of  the  first  kind  [29,31].  Similarly 
the  phase  structure  function  is 


A/,  (d)  = 


3.089  [  — 
r0 


poo 

(  nd\ 

X  / 

1  —  Jo  [  T —  ) 

Jo 

VWJ 

2? xLl  .  (X Lk2 
1  +  A LA  sm 


(A  +  47T2)11/6 


.  (96) 


The  structure  functions  are  plotted  together  in  Fig.  6.  Also  shown  are  the  separa¬ 
tions  required  for  phase  anisoplanatism  (uncorrelated  phase)  and  tilt  anisoplanatism 
(uncorrelated  tilt)  which  are 


00  =  0.949  (k2ClLm)  3/5 , 


(97) 


and 


_  0.319AZT/6 

°TA  ~  CnL3/2 


(98) 


respectively  [18,31].  This  shows  that  two  points  separated  such  that  they  have  log- 
amplitude  independence  still  can  tilt  globally  in  a  similar  manner.  This  allows  a  single 
tracking  mechanism  to  be  used  for  multiple  separated  beams  without  sacrificing  too 
much  wander  between  beams.  The  separation  chosen  for  further  analysis  is  twice 
the  Fresnel  zone  size  of  pc  =  ( L/k )1//2  where  k  =  2ir/\  [1],  For  this  scenario,  the 
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Figure  6.  The  normalized  amplitude  and  phase  structure  functions  are  shown.  When 
the  structure  function  no  longer  increases  with  increased  separation  d ,  then  points 
separated  by  that  distance  are  uncorrelated  from  each  other.  The  vertical  lines  show 
the  separations  due  to  the  isoplanatic  angle  9g,  the  tilt  isoplanatic  angle  Ota-,  and  the 
chosen  separation  of  twice  the  Fresnel  zone  size  2 pc. 


separation  is  2 pc  =  31.4cm. 

2.3.5  Sums  of  gamma-gamma  random  variables. 

The  received  irradiance  from  a  Gaussian  beam  can  be  modeled  as  a  Gamma- 
Gamma  random  variable.  A  multiple-beam  model  can  then  be  treated  as  a  sum  of 
independent  Gaussian  beams,  which  can  be  modeled  as  a  sum  of  independent  Gamma- 
Gamma  random  variables.  The  development  of  the  theoretical  approximation  for 
the  sum  of  multiple,  independent,  identically  distributed  Gamma-Gamma  random 
variables  was  recently  developed  by  Chatzidiamantis  et  al.  It  begins  with  the  sum  of 
the  product  of  two  independent  Gamma  random  variables  [8] 

N 

IN  =  ^2lxiyi,  (99) 

i—  1 
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where  N  is  the  number  of  independent  beams.  This  can  be  rewritten  as 


.  i=l  J  \  i=  1 

JV- 1  M 


(100) 


%—  1  j=i-\- 1 


In  this  form,  the  equation  is  seen  to  be  simply  the  scaled  product  of  the  sum  of  two 
Gamma  random  variables  plus  an  error  term  e  given  by 


N- 1  N 


(101) 


i= 1  j=i- (-1 


Since  the  Gamma  distribution  exhibits  infinite  divisibility  each  sum  of  independent, 
identically  distributed  Gamma  random  variables  is  equivalent  to  a  Gamma  random 
variable  with  shaping  parameters  of  =  Not  and  (3n  —  N/3,  respectively.  The 
two  resulting  Gamma  random  variables  then  have  variances  of  ax/N  and  cry/N.  Us¬ 
ing  Eq.  (57)  and  replacing  a  and  /3  with  cx^r  and  (3n  produces  the  new  multi-beam 
PDF  approximation.  Chatzidiamantis  et  al.  further  refined  the  approximation  by 
accounting  for  the  error  term  by  adding  a  correction  factor  to  ctjv  such  that 


(102) 


CKjv  =  Not  +  e^v- 


The  correction  factor  e n  can  be  approximated  by  minimizing  the  difference  between 
the  first  four  moments  generated  using  the  ejy  with  a  single  Gamma-Gamma  random 
variable  and  the  moments  generated  without  the  correction  factor  using  multinomial 
expansion  for  N  Gamma-Gamma  random  variables.  Using  non-linear  regression,  the 
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result  can  be  closely  approximated  to  be 


ejv  =  (1  -  N) 


0.127  + 0.95a +  0.0058/3 
1  +  0.00124a  +  0.98/3 


(103) 


2.4  Wave-optics  simulations 


Chapter  IV  develops  a  model  for  multiple,  spatially  separated  beams  based  on 
the  work  shown  in  Sections  2.3.4  and  2.3.5.  The  new  model  is  validated  through 
comparison  with  digital  wave-optics  simulation  of  the  beams  propagating  through 
optical  turbulence. 

Wave-optics  simulation  starts  with  the  Fresnel  diffraction  integral  in  Eq.  (19). 
The  diffraction  integral  can  be  converted  into 


U  (x!,yi)  =  U  ( x0,y0 ) 


Akz  ,  ,  ' 

iXz 


(104) 


where  U  (xo,yo)  is  the  input  held  and  x  denotes  the  convolution  of  the  input  held 
with  the  free-space  amplitude  spread  function  [50]. 

At  this  point  Schmidt  introduces  operator  notation  which  is  useful  for  writing  the 
equations  more  compactly  [50].  The  operators  used  are: 


Q[c,  r]  {U  (r)} 
V[b,r]{U{r)} 
F[r,f]{U(r)} 
r-'l f,r]{U(f )} 
n[d,r1,r2]{U  (ri)} 
Q2  [d,r]  {U  (r)} 


eik^2U(r) 


bU  ( br ) 

/OO 

U  (r)  e~i27rf  rd r 

-CO 

/CO 

U  (/)  ei2nf  rd r 

-CO 

i  r°° 

——  /  U  (ri)  el^r2~ri^  dri 

iXa  1 


'  —CO 
^2 2di  |2 


e-2f  HV(r)  =  <2 


47T2 

—d-r 


lU(r)}. 


(105) 

(106) 

(107) 

(108) 

(109) 

(110) 
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Using  this  notation  Eq.  (104)  becomes 


U  (r2)  =  T-'  [r2,  /  J  H  rj  {U  (r,)} ,  (111) 


where  H  (/)  is  the  transfer  function  of  free-space  propagation  given  by 


H  (f  )  =  eikze~i',rXz{fxi+fYl)  ^ 


(112) 


This  technique,  known  as  the  angular-spectrum  method,  is  used  by  the  wave-optics 
simulation  software  WaveProp  (Matlab  toolbox  provided  by  the  Optical  Sciences 
Company)  [57].  To  accommodate  different  grid  resolutions  between  the  input  and 
output  planes,  the  scaling  factor  m  =  5i/52  is  used,  where  and  52  are  the  grid 
spacings  used  in  the  numerical  representation  of  the  input  and  output  planes  respec¬ 
tively.  Using  the  operator  notation  above  and  noting  that  ±m  may  be  used  in  the 
angular-spectrum  form  of  diffraction,  Eq.  (Ill)  becomes 


u  (r2) 
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(114) 


The  simulation  propagates  the  fields  from  one  plane  to  the  next  via  Eq.  (113).  To 
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satisfy  Nyquist  sampling  requirements  the  inequalities  [13,50] 


Sn  < 

N  > 

Azi  < 


Xz  Dn6[ 

Di 

Di  Dn  XAz 

Wl+Wn  +  2^fn 

min  (<5i,  Sn)2  N 
A 


(115) 

(116) 
(117) 


are  used.  The  parameters  <5i  and  5n  are  the  grid  spacing  sizes  at  the  initial  plane 
and  the  final  plane  with  n  —  1  partial  propagations,  Di  and  Dn  are  the  diameters  of 
the  areas  of  interest  at  the  first  and  last  plane,  and  N  is  the  number  of  grid  points 
across  the  field.  To  determine  the  appropriate  propagation  grid  parameter  first  N, 
Si,  and  Sn  are  chosen  based  on  the  first  two  equations  then  choose  the  minimum 
number  of  partial  propagations  by  using  the  third  equation  resulting  in  at  least  n  = 
ceil  (Az/A;smax)  +  1  partial  propagations.  It  should  be  noted  that  the  number  of 
partial  propagations  can  always  be  increased  by  using  smaller  partial  propagation 
distances. 

To  simulate  the  effects  of  turbulence,  random  phase-screens  are  added  to  the 
fields  after  each  propagation  step.  This  collapses  the  phase  effects  of  each  partial 
propagation  into  a  discrete  representation.  To  ensure  that  the  receive  area  of  interest 
captures  the  necessary  amount  of  light,  D\  and  D2  are  replaced  with 


D[ 

XAz 

(118) 

=  Di  +  c 

^0  ,rev 

XAz 

d'2 

=  Do  +  c - , 

r0 

(119) 

where  rojrev  is  the  coherence  diameter  computed  for  light  propagating  in  the  reverse 
direction  and  c  is  an  adjustable  parameter.  Values  of  c  range  from  2  to  8  with  c  =  2 
capturing  97%  of  the  light.  The  sampling  analysis  then  can  be  redone  using  D[  and 
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D'2  in  place  of  Di  and  D2  respectively. 
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III.  Previous  Work 


Research  for  improving  the  effectiveness  of  free-space  laser  communications  sys¬ 
tems  encompasses  several  fields.  The  research  focuses  on  those  fields  most  applicable 
to  an  airborne  FSOC  system.  The  airborne  application  requires  simplicity  and  a  small 
footprint,  both  in  size  and  weight,  while  providing  the  most  availability  over  a  wide 
range  of  conditions.  To  take  full  advantage  of  the  benefits  of  FSOC,  an  integrated 
approach  is  required.  In  this  manner,  each  component  of  current  FSOC  research  is 
surveyed  with  the  intention  of  identifying  areas  of  improvement  or  applicability  to 
the  research. 

The  first  area  of  research  involves  optimizing  the  transmitter  and  receiver  for  the 
the  FSOC  application  when  considering  the  effects  of  atmospheric  turbulence.  Next, 
various  diversity  techniques  are  considered  that  attempt  to  average  out  the  spatial 
and  temporal  variations  at  the  receiver.  After  that,  some  wavefront  control  techniques 
including  lower-order  techniques  such  as  tilt  correction  are  considered.  These  three 
areas  deal  primarily  with  mitigating  the  occurrence  of  intensity  fades  at  the  receiver 
caused  by  atmospheric  turbulence.  In  an  attempt  to  mitigate  the  impact  of  intensity 
fades,  the  use  of  various  modulation  schemes  is  investigated. 

3.1  Transmitter /receiver  design 

The  design  of  the  transmitter  and  receiver  plays  a  large  role  in  the  overall  ef¬ 
fectiveness  of  a  FSOC  system.  The  trade  off  between  between  beam  size  and  beam 
wander  was  shown  by  Churnside  [10].  His  results  showed  that  in  weak  turbulence, 
the  intensity  fluctuations  can  be  reduced  by  adjusting  the  beam  size.  Yenice  et  al. 
proposed  a  scheme  that  adaptively  adjusts  the  beam  size  based  the  observed  tur¬ 
bulence  conditions  [66,67].  By  adjusting  the  beam  size  based  on  the  atmospheric 
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conditions,  the  effects  of  turbulence  can  be  minimized.  More  recently,  it  was  shown 
that  a  linear  array  of  Gaussian  beams  propagated  through  turbulence,  with  spreads 
adjusted  sufficiently,  can  have  the  same  angular  spread  as  a  single  Gaussian  beam 
propagated  through  a  vacuum,  with  the  same  angular  spread,  along  the  dimension 
of  the  array  [25,28].  While  optimizing  the  angular  spread  to  minimize  the  effects 
of  turbulence  is  desirable,  those  benefits  must  always  be  weighed  against  the  cost 
of  increased  beam  size  which  would  decrease  the  average  power  incident  on  the  re¬ 
ceiver  [58]. 

Another  way  to  reduce  the  impact  of  atmospheric  turbulence  on  a  propagated 
beam  is  through  decreasing  the  beam’s  spatial  coherence.  In  addition  to  spreading 
out  the  beam,  the  reduced  coherence  also  reduces  the  intensity  fluctuations  caused 
by  constructive  and  destructive  interference  of  the  beam  with  itself  as  it  propagates 
through  the  turbulence  [45,46].  While  the  speckle  reduction  is  beneficial,  the  increased 
beam  spread  can  result  in  less  power  at  the  receiver,  particularly  for  very  long  paths. 
Recent  work  has  examined  the  respective  trade  offs  and  has  attempted  to  optimize 
the  received  intensity  [9]. 

On  the  receiver  side,  the  most  obvious  method  to  increase  to  total  power  received 
is  to  increase  the  aperture  size.  Additionally,  as  the  aperture  size  increases,  the  fluc¬ 
tuations  decrease,  especially  as  the  size  of  the  receive  aperture  increases  beyond  the 
correlation  width  pc  =  ( L/k )_1//2  for  a  collimated  beam.  For  larger  apertures,  the 
beam  experiences  uncorrelated  intensity  variations  due  to  turbulence  over  the  aper¬ 
ture.  Whether  one  assumes  direct  detection  across  the  aperture,  or  the  incident  light 
is  focused  onto  a  detector  or  coupled  onto  a  fiber,  the  variations  are  averaged  out, 
a  phenomenon  commonly  referred  to  as  aperture  averaging.  The  phenomenon  has 
been  studied  extensively  to  determine  the  performance  improvements  and  the  proba¬ 
bility  distribution  of  the  received  intensity  [11,17,43,69-72].  Zhao  et  al.  showed  the 
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relationship  between  receiver  size  and  propagation  distance  and  determined  optimal 
aperture  sizes  based  on  propagation  distance  and  turbulence  strength  [74],  Other 
recent  work  showed  good  agreement  between  the  lognormal  PDF  and  wave-optics 
simulations  under  weak  scintillation  conditions  while  the  gamma-gamma  PDF  agreed 
with  wave-optics  simulations  under  weak  and  moderate  scintillation  [34].  For  the 
research  described  in  Chapters  IV-VI,  the  aperture  size  is  fixed. 

When  the  incident  light  is  coupled  into  a  fiber,  the  signal  typically  needs  to  be 
amplified,  either  optically  or  electrically.  An  erbium-doped  fiber  amplifier  (EDFA), 
which  optically  enhances  the  signal  but  is  susceptible  to  increased  noise  due  to  spon¬ 
taneous  emission,  increases  the  power  at  the  detector  [22].  Additional  investigations 
showed  that  optical-preamplification  of  30dB  was  strong  enough  to  overcome  ther¬ 
mal  noise  while  still  being  below  the  shot  noise  [42].  Another  method  of  amplifying 
the  signal  uses  an  avalanche  photo  diode  (APD),  which  avoids  noise  resulting  from 
spontaneous  emission  but  introduces  avalanche  noise  [24],  The  research  models  the 
combined  noise  effects  as  AWGN,  as  well  as  direct  detection  models  and  detectors 
without  amplification. 

3.2  Diversity  techniques 

Diversity  techniques  take  advantage  of  the  spatial  and  temporal  variations  in 
atmospheric  conditions  and  average  the  result.  This  includes  using  multiple  trans¬ 
mitters,  multiple  receivers,  and  time  interleaving,  with  the  goal  of  propagating  beams 
through  independent  paths  to  average  out  the  effect  of  turbulence.  Spatial  diversity, 
as  seen  in  the  previous  section,  can  be  accomplished  through  manipulation  of  the 
transmitted  beam  size  or  the  receive  aperture  size  [1,2].  Much  research  has  been 
done  on  multiple  transmitters  and  multiple  receivers  [7,23,35,40,48].  Both  of  these 
techniques  take  advantage  of  the  uncorrelated  spatial  atmospheric  effects  to  average 
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out  the  intensity  variations.  Chatzidiamatis  et  al.  combined  the  use  of  multiple 
transmitters  and  multiple  receivers  as  the  motivation  for  developing  a  model  for  the 
sum  of  N  independent  Gamma-Gamma  random  variables  [8].  Finally,  since  different 
wavelengths  diffract  differently,  over  a  long  path  beams  of  different  wavelengths  can 
be  combined  to  average  out  the  effects  of  turbulence  [39]. 

Recent  work  by  Louthain  and  Schmidt  developed  a  theoretical  basis  for  the  sepa¬ 
ration  required  for  two  points  to  scintillate  independently  [31,33].  Their  work  showed 
that  a  separation  distance  of  one  to  three  times  the  Fresnel  zone  size  pc  would  result 
in  anisoplanatic  scintillation  of  two  points.  The  variation  depends  primarily  on  the 
initial  beam  conditions,  with  a  collimated  beam  having  an  anisoplanatic  separation  of 
2 pc.  The  research  in  the  later  chapters  uses  this  separation  and  applies  it  more  than 
two  transmitters  to  ensure  that  each  transmitter  is  anisoplanatic  from  every  other 
transmitter. 

3.3  Wavefront  control 

Wavefront  control  systems  use  measurements  of  the  incident  light  to  correct  for 
wavefront  aberrations  in  real  time  to  compensate  for  distortions  that  arise  from  prop¬ 
agation  through  the  turbulent  atmosphere.  A  wavefront  control  system  can  be  im¬ 
plemented  on  either  or  both  ends  of  the  transmission  path.  On  the  transmitter  side 
a  reference  beam,  typically  a  point  source,  located  at  the  receiver  end  is  used  to 
measure  the  phase  aberrations  via  a  wavefront  sensor  on  the  transmitter  end.  The 
conjugate  of  the  phase  is  then  imparted  on  the  outgoing  beam  whereby  propagation 
in  the  atmosphere  subsequently  removes  it.  To  implement  corrections  on  the  receiver 
end,  a  reference  at  the  transmitter  is  used,  sometimes  the  transmit  beam  itself.  In 
this  case,  the  conjugate  phase  is  applied  prior  to  the  detector  or  fiber  coupler.  Typ¬ 
ically  wavefront  correction  systems  separate  measurement  and  compensation  for  tilt 
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from  the  higher-order  aberrations.  Tilt  causes  the  beam  to  wander  and  accounts  for 
84%  of  the  total  beam  distortion,  while  the  higher-order  aberrations  cause  the  beam 
to  break  up  [47,49,60]. 

Compensation  for  tilt  aberrations  includes  tracking,  pointing,  and  jitter  control. 
Tracking  refers  to  the  real-time  estimation  on  compensation  to  ensure  that  the  re¬ 
ceived  beam  is  centered  on  the  receiver  and  occurs  on  the  receiver  end.  Pointing 
refers  to  directing  the  beam  in  a  manner  that  ensures  the  outgoing  beam  is  centered 
on  the  target.  Jitter  control  refers  to  compensating  for  the  platform  vibrations  and 
can  occur  at  either  end.  However,  due  to  the  transmission  length,  small  vibrations 
at  the  transmitter  could  result  in  large  deviations  at  the  receiver.  For  example,  a  10 
/irad  deviation  at  the  transmitter  could  result  in  a  1  meter  beam  movement  for  a  100 
km  path.  Recent  research  has  shown  great  capability  with  only  tilt  compensation 
for  ground  to  space  optical  communications  [27].  Crabtree  showed  that  using  a  peak 
tracker,  as  opposed  to  a  traditional  centroid  tracker  or  adaptive  optics,  can  signifi¬ 
cantly  reduce  fade  probability  [14].  Other  techniques  to  compensate  for  jitter  include 
optimizing  the  beam  divergence  and  using  multiple  receive  apertures  with  weighted 
gains  to  minimize  the  BER  [4, 58] .  Most  of  these  techniques  concentrate  on  ground  to 
space  optical  communications  where  tilt  disturbances,  mainly  due  to  jitter,  dominate. 

Higher-order  adaptive  optics  has  been  shown  to  reduce  the  BER  experimen¬ 
tally  [61,62].  The  goal  here  is  to  reduce  scintillation  and  therefore  signal  fades  by 
compensating  for  atmospheric  turbulence.  Additionally,  adaptive  optics  can  be  used 
to  maximize  fiber  coupling  efficiency  [64],  The  research  in  Chapters  IV-VI  attempts 
to  minimize  the  effect  of  atmospheric  turbulence  through  signal  diversity  without  the 
use  of  complicated  higher-order  wavefront  control  schemes. 
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3.4  Modulation  techniques 


The  performance  of  a  FSOC  system  is  not  only  dependent  on  the  power  received, 
but  also  on  the  type  of  modulation  scheme  used.  The  two  broad  categories  of  mod¬ 
ulation  are  coherent  and  incoherent.  Coherent  techniques  encode  the  information  of 
the  light  similar  to  RF  systems  and  require  very  precise  timing.  Incoherent  tech¬ 
niques  vary  the  transmitted  power  without  the  need  of  precise  timing.  Both  types  of 
modulations  have  their  strengths  and  weaknesses. 

The  need  for  precise  timing  in  coherent  techniques  arises  from  the  high  frequency 
of  light.  However,  this  provides  the  potential  for  extraordinary  theoretical  data  rates. 
Lange  et  al.  demonstrated  a  5.625  Gbit/s  link  over  a  142km  terrestrial  path  using 
binary  phase-shift  keying  (BPSK)  [26].  They  used  a  beacon  and  tracking  system 
resulting  in  a  BER  between  ICR4  and  10-6  and  showed  the  potential  robustness  of  a 
coherent  link,  even  in  the  atmosphere.  Other  research  has  shown  that  using  multiple 
transmitters  in  concert  with  a  heterodyne  modulations  scheme  can  also  improve  BER 
performance  [22]. 

Incoherent,  or  direct  detection,  modulates  the  transmitted  power,  the  simplest 
being  OOK  [41].  In  this  binary  digital  communication  scheme,  a  ‘1’  and  ‘O’  are 
represented  by  a  pulse  of  light  or  the  absence  of  a  pulse.  The  advantage  in  simplic¬ 
ity  is  countered  by  the  fact  that  fades  due  to  turbulence  could  result  in  pulses  not 
being  detected.  Variable  threshold  techniques  have  been  proposed  to  address  this 
problem  [33].  Another  technique  is  polarization  shift  keying  (PolSK)  [59,73].  This 
method  takes  advantage  of  the  isotropic  nature  of  turbulence,  which  does  not  affect 
the  polarization  of  the  transmitted  signal.  In  this  case,  the  intensity  fades  due  to 
turbulence  affect  all  symbols  equally.  For  a  binary  signaling  scheme,  this  allows  for 
a  comparison  between  two  signals,  as  opposed  to  a  comparison  against  a  threshold. 
Similarly,  pulse  position  modulation  (PPM)  uses  a  pulse  at  different  temporal  posi- 
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tions  within  a  symbol  and  compares  the  power  detected  at  each  position  [35].  PPM 
takes  advantage  of  the  fact  that  transmitted  symbols  are  much  shorter  in  duration 
than  the  rate  of  change  of  the  turbulence.  The  disadvantage  is  the  more  complex 
timing  and  detection  scheme  required. 

For  those  situations  where  signal  fades  can  completely  disrupt  an  optical  link, 
hybrid  solutions  have  been  investigated  [15].  This  uses  an  optical  link  when  available 
and  switches  to  an  RF  method  when  the  BER  gets  too  high.  For  the  research  de¬ 
scribed  later,  OOK  with  a  single  transmitter  is  used  as  a  baseline,  but  is  compared 
with  OOK  and  PPM  combined  with  multiple  transmitters. 

3.5  Post  processing 

Once  the  symbol  is  received  and  the  bit  estimate  is  determined,  any  resulting  error 
would  need  to  be  detected  or  corrected  to  ensure  effective  communications.  Numer¬ 
ous  forward  error  correction  (FEC)  codes,  originally  developed  for  RF  channels,  have 
been  investigated  for  their  use  with  FSOC.  Various  block,  convolutional,  Turbo,  and 
Reed-Solomon  codes  all  showed  improvements  [38,68,76].  Zhao  et  al.  showed  that  for 
relatively  weak  fading  channels,  the  use  of  Reed-Solomon  encoding  can  improve  BER 
performance  to  nearly  that  of  a  fade-free  channel  [75].  Other  research  investigated  the 
use  of  FEC  with  the  use  of  multiple  transmitters  [3,22],  Because  most  of  these  tech¬ 
niques  were  originally  developed  for  RF  systems  which  primarily  experience  AWGN, 
they  may  not  provide  the  best  solution  for  FSOC.  FSOC  systems  experience  fades 
which  produces  periods  of  reduced  performance.  Any  FEC  or  other  post  processing 
technique  should  take  this  into  account,  minimizing  the  impact  of  fades. 

This  research  described  later  focuses  primarily  on  the  atmospheric  transmission 
and  symbol  detection  aspects  of  the  FSOC  system.  Improvements  in  these  areas 
result  in  a  reduced  BER.  While  the  use  of  FEC  is  not  part  of  this  research,  it  is 
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conceded  that  any  FEC  implementation  on  top  of  the  techniques  considered  would 
only  serve  to  further  reduce  the  overall  BER. 
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IV.  Irradiance  Modelling 


The  effects  of  atmospheric  turbulence,  as  they  pertain  to  a  FSOC  system,  include 
the  effects  of  scintillation  of  the  propagating  beam.  These  intensity  fluctuations  vary 
spatially,  as  well  as  temporally  along  the  propagation  path.  If  the  overall  intensity  of 
the  beam  were  simply  to  increase,  it  would  still  not  address  the  issues  associated  with 
the  deep  fades.  These  deep  fades  ultimately  impact  the  receiver’s  ability  to  estimate 
the  digital  signal,  degrading  overall  communications  system  performance. 

The  first  approach  investigated  here  to  address  the  problem  of  signal  fades  uses 
multiple,  incoherent,  independent  beams  to  “average  out”  the  fades.  The  basic 
premise  is  to  take  advantage  of  multiple  beams,  traveling  through  independent  atmo¬ 
spheric  channels,  which  would  be  unlikely  to  all  experience  fades  at  the  same  time. 
Using  independent  beams  provides  a  best-case  scenario  which  bounds  the  performance 
of  any  multiple  transmitter  system  [55].  This  chapter  begins  with  an  analytical  model 
of  multiple  independent  beams,  and  then  compares  the  model  to  a  wave-optics  simula¬ 
tion  for  a  given  scenario.  Once  the  independent  beam  case  is  developed  the  multiple, 
spatially  separated  beam  (and  correspondingly  partially  correlated)  case  is  addressed. 
First  a  model  for  the  PDF  of  multiple  beams  is  shown,  followed  by  a  comparison  to 
wave-optics  simulations. 

4.1  Independent  beams 

Previously,  the  statistics  of  a  single  aperture-averaged  Gaussian  beam  propagat¬ 
ing  through  turbulence  was  developed.  The  case  of  multiple  independent  provides  a 
limiting  case  for  multiple  beams.  By  using  the  solution  for  multiple,  independent, 
identically  distributed  Gamma-Gamma  random  variables  developed  by  Chatzidia- 
mantis  combined  with  the  shaping  parameters  derived  from  the  aperture  averaged 
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100  km 

Figure  7.  Basic  scenario  description. 


irradiance  shown  in  Eq.  (64)  and  Eq.  (66),  a  bounding  model  for  multi-beam  FSOC 
is  derived.  This  model  bounds  the  performance  for  multiple  transmitters  by  produc¬ 
ing  a  best-case  scenario  where  each  beam  is  completely  independent,  and  therefore 
maximizing  the  averaging  effect  between  beams. 

4.1.1  Simulation  description. 

A  wave-optics  simulation  of  a  nominal  FSOC  system  was  performed  using  the 
angular-spectrum  method  to  evaluate  the  Fresnel  diffraction  integral  [50].  The  simu¬ 
lation  was  designed  to  simulate  uncorrelated  beams  propagating  through  turbulence, 
but  done  in  such  a  way  that  there  was  an  easy  extension  to  partially  correlated 
beams.  The  scenario  and  propagation  grids  were  designed  with  this  ultimate  intent 
in  account. 

The  scenario  for  the  simulation  emulates  a  nominal  FSOC  system  shown  in  Fig.  7. 
A  basic  list  of  the  parameters  used  in  the  simulation  is  shown  in  Table  1.  The 
turbulence  strength  of  =  1  x  10-17m _2/3  equates  to  an  altitude  of  13. 5km  using 
the  Hufnagel-Valley  turbulence  profile. 

4.1.2  Propagation  setup  and  verification. 

The  sampling  requirements  were  determined  using  the  methods  outlined  by  Coy 
to  avoid  aliasing  while  still  adequately  sampling  the  field  in  the  telescope  pupil  [13]. 
Using  Eqs.  (1 15)— (119)  a  Matlab  routine  was  developed  to  simplify  the  sampling 
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Table  1.  Simulation  Parameters. 


Parameter 

Value 

Transmitter  Power,  P 

1W 

Transmitter  Gaussian  Beam  Waist,  Wo 

5cm 

Detector  Diameter,  ddet 

10cm 

Wavelength,  A 

1550nm 

Index  of  Refraction  Structure  Constant,  C2 

1.0  x  10-17m-2/3 

Path  Length,  L 

100km 

analysis  process.  Fig.  8  shows  the  implemented  routine  which  gives  a  contour  plot 
of  the  required  grid  size  for  the  transmitter  and  receiver  grid  spacings.  The  process 
resulted  in  a  chosen  transmitter  grid  spacing  of  24.3mm  and  a  receiver  grid  spacing 
of  7.3mm  along  with  a  required  grid  size  of  29  =  512  points  across  the  grid.  Some 
preliminary  testing  showed  that  using  a  512  x  512  grid  still  introduced  small  amounts 
of  aliasing  so  a  1024  x  1024  grid  (the  next  size  up)  was  used  instead.  Using  the  larger 
grid  size  and  a  propagation  distance  of  100km,  Eq.  (117)  resulted  in  a  minimum  of  4 
propagation  steps. 

To  verify  the  simulation,  a  pair  of  mutually  incoherent  Gaussian  sources  were 
vacuum-propagated  to  the  receive  plane  and  the  resulting  irradiance  was  compared 
with  the  propagated  irradiance  using  Eqs.  (24)-(27).  The  beams  were  separated  by  a 
distance  of  2 pc  =  31. 4cm.  Fig.  9  shows  that  the  simulated  and  analytical  irradiance 
profiles  compare  favorably.  Due  to  the  propagation  distance,  beam  divergence,  and 
small  separation,  the  two  beams  cannot  be  resolved  but  instead  almost  appear  like  a 
single  Gaussian  beam. 

Finally,  the  phase  screen  generation  was  verified  by  comparing  the  phase  structure 
function  of  the  phase  screens  with  the  theoretical  structure  function  of 

/  \  5/3 

D^r)  =6.88  f  y  \  .  (120) 

Since  the  simulation  scenario  used  a  constant  turbulence  strength  along  the  path 
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Figure  8.  Sampling  analysis  for  wave-optics  simulations.  The  figure  is  a  contour  plot 
of  the  required  grid  size  given  the  transmitter  spacing  (along  the  x-axis)  and  receiver 
spacing  (along  the  y-axis).  The  contours  are  shown  as  log2  (IV).  Chosen  grid  spacing 
shows  2.4cm  for  the  transmitter  plane  and  7.3mm  for  the  receiver  plane.  The  analysis 
determined  that  a  grid  size  of  512  x  512  was  sufficient,  however  1024  x  1024  was  ultimately 
used. 
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Figure  9.  Comparison  of  the  simulated  and  analytical  irradiance  profile  for  a  pair  of 
vacuum-propagated  beams. 


48 


r  [m] 

Figure  10.  Comparison  of  the  phase  structure  function  from  the  randomly  drawn  phase 
screens  and  analytical  model.  The  calculated  line  was  generated  from  100  random 
realizations  of  all  phase  screens  spaced  evenly  along  the  path. 


and  the  phase  screens  were  equally  spaced  between  the  transmitter  and  receiver,  each 
phase  screen  has  the  same  correlation  distance  or  r0  value.  Therefore,  the  comparison 
was  done  using  all  of  the  phase  screens  along  the  path.  Fig.  10  shows  the  verification 
of  the  phase  screen  correlation  properties  via  the  structure  function. 

4. 1.2.1  Processing  the  propagated  beams. 

For  this  simulation,  the  initially  collimated  Gaussian  beam  was  transmitted  through 
random  realizations  of  the  atmosphere.  The  resulting  irradiance  at  the  receive  plane, 
collected  as  a  two-dimensional  grid  of  irradiances,  was  masked  based  on  the  size  of 
the  receive  aperture.  Data  points  outside  the  aperture  were  multiplied  by  zero,  while 
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points  within  the  aperture  were  multiplied  by  one.  The  resulting  irradiance  profile 
was  then  integrated  in  both  dimensions  of  the  grid  to  produce  a  single  total  received 
power  across  the  aperture.  The  total  received  power  then  represented  a  single  data 
point.  Another  random  realization  was  produced,  resulting  in  another  data  point, 
and  so  on  creating  n  single-beam  received  powers. 

To  create  multi-beam  results,  the  received  powers  were  combined  in  groups  of 
N  to  correspond  with  an  TV-beam  transmitter.  For  example,  for  dual-beam  results 
the  first  two  received  powers  were  averaged  for  one  data  point,  the  third  and  fourth 
irradiances  were  averaged  for  the  next  data  point,  and  so  on  until  all  of  the  irradiances 
were  considered.  Using  this  method  creates  n/N  data  points  for  N  beams.  The  data 
points  were  then  sorted  from  lowest  to  highest.  By  assigning  each  data  point  a 
marginal  probability  of  N/n  a  CDF  was  created. 

4.1.3  Comparison  of  model  and  simulation. 

A  comparison  between  the  aperture-averaged  Gamma-Gamma  CDF  from  Eq.  (64) 
using  a  and  f3  resulting  from  Eqs.  (58)  and  (66)  with  a  wave-optics  simulation  for 
a  single  beam,  nominal  FSOC  system  is  shown  in  Fig.  11.  The  single-beam  results 
show  that  the  well-established  analytical  model  for  a  single  beam  produces  a  good 
fit  to  the  simulation  data.  This  further  confirms  that  the  simulation  is  sufficiently 
modeling  the  engagement  scenario. 

The  results  shown  in  Fig.  12  for  one  to  seven  beams  show  good  agreement  with 
the  model  and  also  demonstrate  a  decrease  in  fade  probability  as  the  number  of 
transmitters  are  increased.  It  also  shows  that  the  marginal  improvement  decreased 
as  the  number  of  beams  used  increased.  The  improvement  between  one  beam  and 
two  beams  at  a  cumulative  probability  of  10-3  is  3.02dBm.  However,  the  difference 
between  two  and  four  beams  is  only  1.95dBm  and  the  difference  between  four  and 
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Received  Irradiance  [dBm] 


Figure  11.  Analytical  and  wave-optics  generated  CDF  for  an  aperture-averaged,  single 
beam  irradiance. 

seven  beams  is  only  1.08dBm.  Using  multiple  transmitters  appears  to  be  an  effective 
way  of  reducing  fades,  but  there  seems  to  be  a  point  of  diminishing  returns. 

The  uncorrelated  beam  case  provides  a  useful  best-case  scenario  for  multi-beam 
transmitters.  Exploring  the  effects  of  and  establishing  a  model  for  partially  correlated 
beams  adds  further  realism  to  this  study. 

4.2  Partially  correlated  beams 

Practical  application  of  a  multi-beam  transmitter  requires  some  finite  spatial  sep¬ 
aration  between  the  beams.  This  means  that  the  beams  would  not  be  independent 
from  each  other,  but  partially  correlated.  This  section  develops  another  model  to  ad¬ 
dress  this  partially  correlated  case.  Following  that,  another  comparison  to  simulation 
is  performed. 
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Figure  12.  Analytical  CDF,  solid  lines,  and  wave-optics  results,  dashed  lines.  Results 
are  shown  for  one,  two,  four,  and  seven  beams  based  on  the  number  of  beams  used  in 
the  spatially  separated,  multiple  beam  configuration. 

4.2.1  Beam  configuration. 

To  minimize  the  amount  of  area  occupied  by  the  multiple  transmitters,  a  hexag¬ 
onal  close-pack  spacing  is  used  [54,56].  Fig.  13  shows  the  arrangement  and  relative 
intensities  of  each  of  the  patterns.  Each  pattern  used  is  centered  in  relation  to  the 
single  beam.  The  closest  any  two  beams  can  be  is  the  chosen  separation  distance,  in 
this  case  2 pc.  For  the  larger  number  of  transmitters,  particularly  five  or  more,  the 
maximum  separation  between  the  furthest  spaced  beams  is  twice  the  chosen  sepa¬ 
ration  distance  or  4 pc.  Those  beams  begin  to  wander  independently,  since  they  are 
spaced  further  than  the  spacing  required  for  tilt  isoplanatism,  which  would  limit  the 
ability  to  use  a  single  tracking  sensor  for  beam  steering.  A  simple  solution  would 
be  to  track  each  beam  individually  at  the  cost  of  increased  system  complexity.  Ad¬ 
ditionally,  regardless  of  the  number  of  transmitters  used,  the  average  power  is  held 
constant,  which  in  this  scenario  is  P  =  1W. 
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Figure  13.  The  configuration  used  for  multiple  transmitters.  Each  configuration  is  a 
subset  of  a  hexagonal  close-pack  grid.  The  top  row  shows  one,  two,  three,  and  four 
transmitters  from  left  to  right,  while  the  bottom  row  shows  five,  six,  and  seven  trans¬ 
mitters.  In  all  cases,  the  total  average  power  across  all  transmitters  is  held  constant. 


4.2.2  Model  development. 

To  account  for  the  finite  beam  separation,  a  method  needed  to  be  developed  to 
account  for  the  partially  correlated  scintillation  resulting  from  several  beams.  The 
independent,  identically  distributed  shaping  parameters  can  be  rewritten  as 
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(121) 

(122) 


which  treats  the  modification  to  a  single  beam  shaping  parameter  as  a  simple  gain. 
There  are  then  two  limiting  cases  that  were  considered,  perfectly  correlated  beams 
and  perfectly  uncorrelated  beams.  For  perfectly  correlated  beams  the  shaping  pa¬ 
rameter  for  multiple  beams  should  be  the  same  as  for  a  single  beam.  For  perfectly 
uncorrelated  beams  the  shaping  parameter  for  multiple  beams  should  be  the  same  as 
for  uncorrelated  beams.  By  raising  the  gain  to  a  function  of  the  correlation  coefficient 
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the  partially  correlated  shaping  parameters  become 


Ctj\r  — 


Na  +  e 


a 


f(p) 


ft,  =  /3N'm. 


(123) 

(124) 


To  satisfy  the  limiting  cases  then 


lim/(p)  =  1, 

p— 1-0 


(125) 


and 


lirn  /  (p)  =  0. 
i 


(126) 


The  correlation  coefficient  itself  is  a  function  of  beam  separation  resulting  in 
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and  has  limiting  cases  of 

lim  p(d)  =  1, 

d — ^0 

and 

lim  p(d)  =  0. 

d^-oo 


(129) 


(130) 


The  correlation  coefficient  p  (d)  was  determined  using  wave-optics  simulations  of 
Gaussian  beam  pairs  of  various  separation  distances  propagated  over  the  simulated 
distance  and  integrated  over  a  common  aperture  located  between  the  propagated 
beams  at  the  receive  plane.  The  resulting  relationship  between  separation  distance 
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and  correlation  coefficient  was  determined  via 


p(d) 


{{I A  (-d/2,  0)  -  jilA  (-d/2,  0))  (IA  (d/2, 0)  -  jilA  (d/2,  0))) 
alA(-d/2,0)alA(d/2,0) 


(131) 


where  (•)  is  the  expectation  operator  and  I  a  (x,  y)  is  the  received  irradiance  inte¬ 
grated  over  the  receive  aperture  for  a  Gaussian  beam  centered  on  x,  y  and  /i/A  (x,  y ) 
and  aiA  (x,  y)  are  the  mean  and  variance  of  I  a  (x,  y)  respectively.  The  integrated 
irradiance  for  each  beam  was  determined  using  Eq.  (90)  with  the  short-term  beam 
waist  Wst  replaced  by  the  long-term  beam  waist  Wlt-  The  simulated  results  are 
shown  in  Fig.  14  along  with  first-order  results  for  a  single-parameter  functional  fit  of 


p  (d)  ~  exp  (— 0.6875d/pc) . 


(132) 


The  value  of  —0.6875  in  the  model  was  determined  by  minimizing  the  integrated  mean 
squared  difference  between  the  model  and  wave-optics  results.  This  functional  form 
also  provides  the  appropriate  limiting  conditions  as  d  approaches  zero  and  infinity. 

The  same  simulation  data  were  then  used  to  produce  CDFs  of  the  integrated 
irradiance  over  the  receive  aperture  for  two  beams  at  various  separations.  Analytical 
CDFs  using  the  shaping  parameters  in  Eqs.  (127)  and  (128)  were  then  produced  using 
an  value  for  /  (p  (d))  which  minimized  in  integrated  mean  squared  difference  between 
the  analytical  CDF  and  the  wave-optics  produced  CDF.  The  resulting  f(p(d))  is 
plotted  against  separation  distance  between  the  beams  in  Fig.  15.  A  model  of  (p(d)) 
is  also  shown  in  Fig.  15  such  that 


(133) 


where  the  value  of  1.4894  was  determined  by  minimizing  the  integrated  mean  squared 
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d/ pc 

Figure  14.  Correlation  coefficient  versus  beam  separation  for  two  Gaussian  beams 
propagated  over  100km  path  with  a\  =  1.0. 

difference  between  the  model  and  simulation  results.  Combining  Eqs.  (132)  and  (133) 
results  in 

f  (p  (d))  «  [1  -  exp  (-0.6875d/pc)]14894  .  (134) 

The  simulation  results  in  Fig.  15  show  that  at  a  separation  distance  greater  than 
3 pc  the  value  of  f(p(d))  begins  to  decrease.  The  aperture-averaged  model  for  the 
Gamma-Gamma  shaping  parameter  assumes  the  receive  aperture  is  centered  on  the 
beam  center.  As  the  beam  separation  increases,  this  assumption  becomes  less  valid. 
To  investigate  the  impact  of  the  displacement  on  the  CDF,  several  normalized  CDFs 
were  created  from  simulation  data  with  varying  receiver  displacements.  The  resulting 
plot  shown  in  Fig.  16  indicates  that  there  is  a  slight  ~  O.ldBm  shift  in  the  CDF 
at  a  probability  of  10”2.  This  error  source  is  small  when  compared  to  the  CDF 
shift  associated  with  increasing  the  number  of  transmitters,  which  is  on  the  order 
of  l.OdBm.  The  results  show  that  the  Gamma-Gamma  CDF  used  for  single-beam 
propagation  adequately  models  off  axis  irradiance  statistics  as  well. 
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Figure  15.  Exponent  value  /  (p(d))  in  Eqs.  (127)  and  (128)  plotted  as  a  function  of  beam 
separation  for  two  Gaussian  beams  propagated  over  100km  path  with  af  =  1.0. 
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Figure  16.  Normalized  CDFs  for  integrated  irradiance  with  receivers  placed  at  various 
off-axis  distances. 
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Figure  17.  Mean  integrated  irradiance  plotted  as  a  function  of  receiver  distance  from 
beam  center. 


Another  potential  source  of  the  decreasing  value  of  /  (p)  beyond  3 pc  in  Fig.  15 
is  the  variation  in  mean  received  irradiance  as  a  function  of  receiver  displacement 
from  beam  center.  Fig.  17  demonstrates  that  as  the  receiver  displacement  increases, 
the  mean  irradiance  decreases.  This,  coupled  with  a  corresponding  decrease  in  the 
variance  as  displacement  is  increased,  provides  a  partial  explanation  for  the  downward 
trend  evident  in  Fig.  15.  However,  the  primary  region  of  interest  lies  between  lpc 
and  3 pc  where  there  is  sufficient  beam  separation  to  average  the  irradiance  between 
beams  but  still  close  enough  for  practical  applications.  At  this  separation  distance 
there  is  still  reasonable  agreement  between  the  simulation  data  and  the  model  shown 
in  Fig.  15. 

Finally,  the  model  also  assumes  that  each  beam  has  the  same  separation  distance 
from  every  other  beam.  While  this  is  true  for  two  or  three  beams,  this  assumption  is 
not  true  for  four  or  more  beams.  However,  the  beam  configuration  shown  in  Fig.  13 
ensures  that  beams  are  no  closer  than  the  chosen  separation  distance.  Additionally 
the  correlation  coefficient  as  a  function  of  separation  distance  shown  in  Fig.  14  indi- 
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X  position  [m] 


Figure  18.  The  union  of  the  seven  beam  configurations  used  for  multiple  transmitters. 

cate  that  correlation  decreases  as  separation  distance  increases.  This  indicates  that 
the  model  in  Eqs.  (127)  and  (128)  produce  a  pessimistic  estimate  for  multi-beam 
performance. 

4.2.3  Simulation  approach. 

The  simulation  configuration  used  for  the  independent  beam  case  was  also  used 
for  the  multiple,  partially  correlated  beam  case.  To  propagate  the  multi-beam  con¬ 
figurations,  a  matrix  of  seventeen  unique  beam  positions  was  created  by  determining 
the  union  of  the  seven  beam  configurations  shown  in  Fig.  18.  After  propagating  the 
beams  and  integrating  them  across  the  receiver,  an  L  x  M  matrix  of  irradiances  was 
created  where  L  is  the  number  of  beam  positions  and  M  is  the  number  of  random 
realizations.  An  M  x  N  evaluation  matrix  was  created  to  convert  the  beam  posi¬ 
tion  irradiances  to  transmitter  configuration  irradiances  where  N  is  the  number  of 
transmitters. 

One  issue  associated  with  the  placement  of  multiple  beams  resulting  from  the 
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necessarily  non-uniform  separation  distance  for  four  or  more  beams  is  the  variation 
in  relative  intensity  over  the  receive  aperture  from  each  of  the  transmitters.  When  the 
beams  are  transmitted  parallel  to  each  other,  beams  placed  further  from  the  center 
clearly  result  in  a  smaller  contribution  to  the  aperture-averaged  irradiance  since  the 
aperture  intercepts  the  beam  edge  and  not  its  center.  As  a  result,  fades  from  beams 
near  the  center  could  have  a  greater  impact  on  the  frequency  of  fades  than  those 
farther  away.  The  model  for  partially  correlated  beams  also  assumes  that  each  beam 
contributes  to  the  mean  irradiance  equally.  To  compensate,  the  transmitted  power 
for  each  beam  was  adjusted  to  ensure  the  mean  power  received  from  each  beam  is 
the  same.  Each  transmitter’s  mean  contribution  to  the  integrated  irradiance  was 
determined  by  using  Eq.  (90)  with  the  short  term  beam  waist  Wst  replaced  by  the 
long  term  beam  waist  Wlt ■  This  resulted  in  an  evaluation  matrix  of 
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Without  compensating  for  the  different  beam  powers  the  evaluation  matrix  would 
simple  contain  values  where  each  column  would  only  be  populated  by  either  a  zero 
or  1/N.  This  effectively  would  allow  the  contribution  from  beam  to  be  equal  for  a 
given  configuration. 

4.2.4  Comparison  of  the  results. 

Using  a  combination  of  multiple  beams  and  adjusting  the  power  among  the  beams 
such  that  the  mean  power  received  from  each  is  the  same,  a  wave-optics  simulation  was 
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done  showing  the  marginal  improvement  as  the  number  of  beams  is  increased.  The 
CDF  results  are  shown  in  Fig.  19.  Additionally,  the  theoretical  estimates  determined 
using  Eqs.  (57),  (127),  and  (133)  were  compared  to  the  CDFs  from  simulation.  The 
reduced  fade  probability  is  apparent  as  the  number  of  transmitters  is  increased.  For 
example,  if  the  fade  threshold  is  determined  to  be  OdBm,  the  probability  of  fade 
for  one,  two,  four,  and  seven  transmitters  is  9.4  x  10~3,  2.0  x  10”3,  3.7  x  10~4,  and 
1.5  x  ICC4,  respectively.  This  shows  that  using  more  than  four  transmitters  does  little 
to  improve  the  overall  performance. 

Fig.  19  shows  that  the  model  provides  a  good  description  of  wave-optics  results  for 
up  to  four  transmitters.  As  the  number  of  transmitters  increases  past  four,  the  model 
for  the  shaping  parameters  in  Eqs.  (127)  and  (128)  appears  to  become  less  valid.  The 
assumption  of  Gamma-Gamma  statistics  for  off-center  beams  is  likely  failing  again. 
Referring  back  to  Fig.  13,  one  can  see  that  as  the  number  of  transmitters  increases, 
separations  between  the  center  of  the  beam  configurations,  where  the  receiver  is  lo¬ 
cated,  and  the  most  distant  beams  increase.  This  may  account  for  the  disagreement 
between  the  model  and  wave-optics  simulations  when  the  number  of  transmitters  is 
greater  than  four. 

Overall,  the  model  provides  a  good  tool  for  estimating  the  fade  probability  for  a 
multiple-transmitter  system.  This  would  give  engineers  a  method  to  investigate  the 
design  trade  offs  associated  with  the  multi-beam  systems. 
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Figure  19.  Analytical  CDF,  solid  lines,  and  wave-optics  results,  dashed  lines,  plotted 
for  one  through  seven  spatially  separated  beams  using  the  configurations  shown  in 
Fig.  13. 
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V.  Communications  System  Performance 


This  chapter  discusses  the  BER  results  for  and  ideal  fixed  threshold  OOK  and 
PPM.  First,  an  analytical  model  for  the  BER  for  a  multi-beam  FSOC  is  developed. 
The  next  section  provides  a  description  of  how  the  irradiance  from  the  wave-optics 
simulation  was  used  in  combination  with  traditional  BER  models  to  calculate  BER. 
Finally,  a  Monte  Carlo  analysis  of  the  BER  was  performed  using  the  irradiances  from 
the  wave-optics  simulation.  All  three  of  these  methods  were  compared  with  each 
other  showing  the  significant  increase  in  performance  when  multiple  beams  are  used 
in  combination  with  PPM  over  baseline  method  of  a  single  beam  with  OOK. 

The  difference  in  performance  between  PPM  and  OOK  highlights  the  difference 
between  using  a  fixed  threshold  OOK  modulation  scheme  which  is  highly  sensitive 
to  deep  signal  fades  versus  a  scheme  which  is  not  significantly  impacted  by  deep 
fades.  In  the  case  of  OOK,  a  deep  fade  due  to  scintillation  causes  the  received  signal 
to  drop  below  the  threshold,  resulting  in  every  received  bit  being  estimated  as  a 
signal  low.  Assuming  the  likelihood  of  a  transmitted  low  or  high  signal  is  equal, 
this  amounts  to  a  50%  BER.  Binary  PPM  places  all  of  the  energy  either  in  the  first 
half  of  the  transmitted  symbol  or  the  second  half.  The  receiver  then  subtracts  the 
received  energy  for  the  first  half  of  the  symbol  from  the  second  half  resulting  in  either 
a  positive  or  negative  value.  With  a  threshold  of  zero,  the  negative  values  correspond 
with  a  signal  low  and  positive  values  with  a  signal  high.  PPM  does  not  require  the 
constant  detection  and  update  of  the  threshold  value,  which  draws  irradiance  from 
the  communications  detector  [33].  Additionally,  PPM  always  has  either  a  rising  or 
falling  edge  in  the  center  of  the  pulse  which  could  be  used  to  aid  in  timing  and 
synchronization  of  the  demodulator  [53].  The  price  of  these  benihts  is  a  doubling 
of  the  bandwidth  necessary  to  transmit  and  receive  the  PPM  symbol  compared  to 
OOK. 
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5.1  Analytical  model  of  BER 


The  analytical  BER  model  starts  with  the  calculation  of  BER  for  a  binary  signal¬ 
ing  scheme.  Using  the  detector  model  Eq.  (1)  and  assuming  that  the  channel  gain  is 
constant  over  the  duration  of  a  symbol,  the  BER  for  a  given  channel  gain  PB  [h]  can 
be  determined  via  [52] 

PB  [h]  =  P  [s0]  P  [Hi\s0l  h}  +  P  [Sl]  P  [H0\Sl,  h] ,  (136) 

where  h  represents  the  channel  gain  over  the  symbol  duration  and  s0,i  represents  the 
two  transmitted  symbols  for  the  binary  signaling  scheme.  P  [so.i]  is  the  probability 
that  s0  or  si  was  transmitted  respectively,  and  P  [//0,ilso,i;  h]  is  the  probability  that 
the  receiver  hypothesized  si  when  so  was  transmitted  and  vice  versa  for  a  given 
channel  gain.  The  assumption  that  h  (t)  is  constant  over  the  symbol  duration  requires 
that  T  <C  1//g  where  T  is  the  symbol  duration  and  fc  is  the  Greenwood  frequency 
from  Eq.  (74).  For  all  practical  systems  of  interest  this  requirement  is  met  [32,33,51]. 
It  should  also  be  noted  that  h  (t)  is  not  a  zero-mean  process  and  is  characterized  using 
the  same  model  developed  in  Chapter  IV  with  a  transmitted  irradiance  of  1W. 

Using  the  model  in  Eq.  (1)  assuming  AWGN,  the  probability  of  the  two  error 
conditions  in  Eq.  (136)  become 
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where  zt  is  the  threshold  from  Eq.  (8),  /z0,i  is  the  test  statistic  mean  values  for  s0  and 
si  respectively,  and  az  is  the  standard  deviation  of  the  test  statistic  for  so  and  si  due 
to  the  AWGN.  The  operator  Q  (•)  is  the  Q-function  defined  as 


OO 


X 


To  calculate  the  total  BER  over  all  possible  channel  gains  for  a  specific  atmospheric 
path,  the  BER  for  each  channel  gain  is  combined  with  the  probability  of  each  channel 
gain  value  via 

OO 

Pb  =  J  PB[h]p(h)  dh,  (140) 

o 

where  p  (h)  is  the  PDF  the  channel  gain.  The  PDF  of  the  channel  gain  equivalent  to 
the  PDF  of  irradiance  in  Eq.  (57)  assuming  a  1W  transmitted  power.  The  integral 
in  Eq.  (140)  was  evaluated  numerically. 

The  difference  between  the  performance  of  an  ideal  fixed-threshold  OOK  modu¬ 
lation  scheme  and  a  PPM  modulation  scheme  simplifies  to  choosing  the  appropriate 
reference  signal  sref  (t).  For  OOK  the  reference  signal  is 


.  T 

Sref  (t)  =  u(t)  0  >t< 


(141) 


while  for  PPM  the  reference  signal  is 


T 

Sref  ( t )  =  — U  ( t )  +2 u(t  —  T / 2)  0  >  t  <  — , 


(142) 


where  u  ( t )  is  the  unit  step  function.  For  the  fixed  threshold  OOK  case  the  threshold 
was  based  on  Eq.  (8).  The  expected  test  statistic  value  for  receiving  an  on  bit,  pi  in 
Eq.  (8),  was  calculated  with  a  channel  gain  determined  from  Eq.  (90)  and  replacing 
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the  short-term  beam  waist  with  the  long-term  waist  found  in  Eq.  (91).  The  variance 
was  calculated  by  summing  the  AWGN  variance  with  the  variance  calculated  from 
Eq.  (59)  using  the  shaping  parameters  from  Eq.  (127). 

For  PPM,  since  /x0  =  —  Hi  the  threshold  is  zero.  This  means  that  regardless  of  the 
channel  gain,  the  threshold  is  always  optimal,  resulting  in  a  better  BER  compared  to 
using  fixed  threshold  OOK.  The  price  for  this  improved  performance  is  an  increase 
in  the  bandwidth  required  by  a  factor  of  two. 

5.2  BER  with  wave-optics 

Calculation  of  the  BER  from  the  wave-optics  results  simply  used  the  channel 
gains  from  each  random  realization  of  the  atmosphere  along  the  propagation  path 
to  calculate  a  BER  based  on  random  draws  of  the  aperture  averaged  irradiace  at 
the  receiver.  The  BER  for  all  the  random  realizations  of  the  channel  gains  were 
averaged  together  to  produce  the  overall  BER  for  the  propagation  path.  The  ideal 
fixed  threshold  was  determined  by  using  the  mean  and  variance  of  the  test  statistic 
from  the  simulation  data  and  changed  depending  on  the  number  of  transmitters  and 
the  AWGN  power.  The  PPM  simulation  used  a  threshold  of  zero.  The  simulation 
results  from  Section  4.2  were  used  to  generate  the  channel  gains  by  calculculating  the 
ratio  between  the  input  irradiance  and  the  integrated  irradiance  across  the  receiver. 

5.3  Monte  Carlo  simulation  of  BER  measurement 

Additionally  a  Monte  Carlo  simulation  was  performed  at  each  time-step  to  validate 
the  analytical  model.  First,  the  Monte  Carlo  simulation  generated  a  random  bit 
pattern  at  an  information  bit  rate  of  lGbps.  The  bit  pattern  was  then  modulated 
at  a  sampling  rate  of  200GHz  using  both  OOK  and  PPM  to  create  the  symbols. 
Following  the  model  in  Eq.  (1),  the  symbols  were  multiplied  by  the  channel  gain 
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from  an  instance  of  the  wave-optics  simulation  and  random  draws  AWGN  noise  were 
added.  Then,  the  modeled  detector  output  was  multiplied  by  the  reference  signal  and 
integrated  over  the  symbol  duration  to  obtain  the  test  statistic  shown  in  Eq.  (2).  The 
test  statistic  was  then  compared  to  the  threshold  to  produce  the  bit  estimates.  The  bit 
estimates  were  then  compared  to  the  original  random  bit  pattern  and  differences  were 
counted  as  bit  errors.  The  process  was  repeated  for  each  of  7,  250  random  realizations 
of  the  wave-optics  simulation.  Prior  to  executing  the  simulation,  the  Monte  Carlo 
simulation  was  validated  using  a  fixed  channel  gain  and  only  AWGN  which  produced 
results  identical  to  the  established  models  [52]. 


5.4  Comparison  of  results 


A  comparison  of  the  results  for  all  three  methods  are  shown  in  Figs.  20  and  21.  The 
calculated  BER  is  shown  as  a  function  of  E^/Nq  where  E b  is  the  transmitted  signal 
energy  per  information  bit  and  N0  is  the  AWGN  power  over  the  receiver  bandwidth. 
Eb/No  is  related  to  SNR  by 


E7  _  SW_ 

Ao  ““  N~R 


(143) 


where  S  is  the  signal  power,  N  is  the  noise  power,  W  is  the  bandwidth  of  the  receiver, 
and  R  is  the  information  bit  rate.  The  results  are  limited  to  the  aperture  size  used  in 
the  simulation,  D  =  10cm,  but  allows  for  a  meaningful  and  fair  comparison  between 
the  various  transmitter  configurations,  since  different  transmitter  configurations  result 
in  different  mean  received  irradiance. 

The  use  of  additional  transmitters  clearly  improves  the  BER  seen  in  Figs.  20 
and  21.  However,  the  improvement  is  more  dramatic  for  the  fixed  threshold  OOK 
than  the  improvement  for  PPM.  By  way  of  comparison  at  a  BER  of  ICC3,  the  gain 
between  one  to  two  transmitters  is  3.4dB  for  fixed  threshold  OOK  and  only  1.6dB 
for  PPM.  Similarly,  the  gain  for  going  from  two  to  four  transmitters  is  2.4dB  for 


Figure  20.  BER  for  multiple  beams  using  a  fixed  threshold  OOK.  The  markers  indicate 
results  from  the  Monte  Carlo  simulation  of  the  communications  link.  The  dashed  lines 
are  the  “analytical”  results  using  the  data  from  the  wave-optics  simulation.  The  solid 
line  is  the  analytical  result  obtained  from  the  newly  derived  model  of  the  received 
intensity.  Results  are  shown  for  one,  two,  four,  and  seven  transmitters. 


Eb/N0  [dB] 

Figure  21.  BER  for  multiple  beams  using  PPM.  The  markers  indicate  results  from 
the  Monte  Carlo  simulation  of  the  communications  link.  The  dashed  lines  are  the 
“analytical”  results  using  the  data  from  the  wave-optics  simulation.  The  solid  line  is 
the  analytical  result  obtained  from  the  newly  derived  model  of  the  received  intensity. 
Results  are  shown  for  one,  two,  four,  and  seven  transmitters. 
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fixed  threshold  OOK  and  only  1.3dB  for  PPM.  While  gain  is  greater  when  fixed 
threshold  OOK  is  used,  the  performance  still  does  not  surpass  that  of  PPM,  which  at 
a  BER  of  1CT3  PPM  is  15.7dB  better  even  when  comparing  seven  transmitters  using 
fixed  threshold  OOK  against  a  single  transmitter  using  PPM.  Figure  20  also  shows 
the  more  dramatic  impact  the  modeling  errors  have  in  estimating  the  BER  for  fixed 
threshold  OOK  as  the  number  of  transmitters  is  increased.  While  the  BER  curve  for 
one  or  two  transmitters  is  modeled  well,  the  assumptions  used  for  determining  the 
PDF  for  additional  transmitters  more  dramatically  impacts  BER  estimates.  This  is 
less  evident  for  the  PPM  case  shown  in  Fig.  21,  where  there  is  much  better  agreement 
between  the  model  and  simulation  results. 

Since  PPM  requires  twice  the  bandwidth  compared  to  OOK,  a  better  comparison 
would  be  to  use  the  extra  bandwidth  to  implement  FEC.  An  Extended  Golay(24,12) 
algorithm  was  used  in  a  Monte  Carlo  simulation  to  generate  24  channel  bits  for  every 
12  data  bits.  The  encoded  bits  were  then  propagated  through  the  channel  at  a  rate  of 
2Gbps,  twice  the  rate  of  the  information  bits.  This  results  in  the  encoded  OOK  signal 
using  the  same  bandwidth  as  PPM.  After  propagation  through  the  channel  and  the 
addition  of  AWGN,  the  channel  bit  estimates  were  decoded  to  produce  information 
bit  estimates.  Errors  between  the  transmitted  information  bits  and  the  information 
bit  estimates  were  counted  as  errors.  Fig.  22  shows  the  results  of  this  simulation. 
For  this  FSOC  scenario,  where  the  bit  errors  for  OOK  are  dominated  by  the  channel 
fades,  fixed  threshold  OOK  with  FEC  shows  little  performance  improvement  over  the 
baseline  of  fixed  threshold  OOK. 
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Eb/N0  [dB] 

Figure  22.  Simulation  results  for  fixed  threshold  OOK  and  PPM  compared  with  FEC 
OOK  The  results  show  that  for  this  scenario,  where  the  channel  gain  of  the  atmosphere 
dominates  the  bit  errors,  FEC  does  little  to  improve  overall  BER.  Results  are  shown 
for  one,  two,  four,  and  seven  transmitters. 
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VI.  Lab  Experimentation 


The  results  from  a  bench-top  experiment  using  an  atmospheric  turbulence  sim¬ 
ulator  (ATS)  were  compared  with  the  model  developed  in  Chapter  IV.  The  in-lab 
experiment  simulated  a  multi-beam  communications  system  by  transmitting  a  colli¬ 
mated  laser  through  the  turbulence  simulator  onto  an  imaging  system.  The  imaging 
system  is  used  to  determine  fade  statistics,  which  were  then  compared  to  the  previ¬ 
ously  developed  model. 

6.1  Lab  description 

To  more  accurately  explore  the  complexities  involved  with  an  actual  FSOC  system, 
a  bench-top  laboratory  demonstration  was  constructed  taking  advantage  of  existing 
equipment.  The  experiment  validates  the  analytical  model  and  wave-optics  simula¬ 
tion  by  measuring  the  aperture-averaged  irradiance  using  multiple  transmitters.  The 
overall  system  configuration  is  described  below,  followed  by  a  more  detailed  descrip¬ 
tion  of  ATS. 

6.1.1  Configuration. 

The  basic  layout  of  the  experimental  system  is  shown  in  Fig.  23.  The  laser  diode, 
a  Thorlabs  LPS-PM1550-FC  mounted  to  a  Thorlabs  TCLDM9  Temperature  Con¬ 
trolled  Laser  Diode  Mount  shown  in  Fig.  24,  emits  a  polarized  1550nm  beam  through 
a  polarization-maintaining  fiber.  The  emitted  power  level  is  tunable  via  a  Thorlabs 
ITC8102  controller  mounted  in  a  PR08000  chassis.  The  polarization-maintaining 
fiber  connects  to  a  Thorlabs  F230FC-1550  fiber  collimator  resulting  in  a  0.9mm  colli¬ 
mated  beam.  The  beam  then  travels  through  the  SAIC  ATS  described  below.  After 
exiting  the  ATS,  the  beam  then  relays  off  of  a  total  of  four  relay  mirrors  used  to 
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IR  Camera 

Figure  23.  Lab  experiment  configuration. 


Figure  24.  (a)  Thorlabs  LPS-1550-FC  Laser  Diode,  (b)  Thorlabs  TCLDM9  Temper¬ 
ature  Controlled  Laser  Diode  Mount,  (c)  Thorlabs  F230FC-1550  1550nm  Fiber  Colli¬ 
mating  Package. 


increase  the  distance  between  the  ATS  and  the  imaging  system.  The  relay  extends 
the  total  path  length  from  a  nominal  lrn  beyond  the  ATS,  used  for  adaptive  optics 
experiments,  to  6.2m  beyond  the  ATS. 

The  imaging  system  consists  of  two  lenses  configured  as  a  4-f  system  where  a  l/4th 
magnification  is  imposed  on  the  beam.  The  resulting  image  is  received  by  a  FL1R 
Alpha  near  infra-red  (NIR)  camera.  The  first  lens  has  a  focal  length  of  300mm,  while 
the  second  lens  has  a  focal  length  of  75mm.  Both  lenses  are  anti-reflection  coated 
for  the  1550nm  beam.  To  find  the  camera  image  plane,  a  picture  was  placed  on  a 
filter  mount  and  the  camera  position  was  adjusted  until  the  picture  came  into  focus. 
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Figure  25.  Imaging  system  showing  the  300mm  lens  and  the  75mm  lens  followed  by 
the  FLIR  Alpha  NIR  camera. 


The  camera  was  locked  down  and  the  picture  removed,  leaving  the  filter  mount  to 
mark  the  camera  imaging  plane  position.  The  whole  system  could  then  be  moved 
together  to  adjust  the  path  length  between  the  imaging  system  and  the  ATS  without 
disturbing  the  imaging  plane  location  relative  to  the  camera. 

Once  the  data  are  captured,  a  digital  mask  is  applied  to  the  saved  images  to 
simulate  various  aperture  sizes.  By  moving  the  collimated  beam,  causing  the  beam  to 
travel  through  a  different  portion  of  the  turbulence  generated  by  the  ATS,  the  system 
can  simulate  multiple  off-axis  beams.  The  images  are  combined  in  post-processing  to 
simulate  multiple  beams  propagating  at  the  same  time. 

6.1.2  Atmospheric  turbulence  simulator. 

The  turbulence  simulator  is  an  SAIC  ATS,  model  0003  shown  in  Fig.  26  [36,44], 
Using  a  combination  of  lenses  and  phase  wheel  modules,  the  ATS  simulates  atmo¬ 
spheric  propagation  for  a  variable  range  of  turbulence  parameters.  By  adjusting  the 
phase  wheel  location,  the  density  of  the  variations  across  the  beam  can  be  changed, 
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Figure  26.  Atmospheric  Turbulence  Simulator. 


since  the  phase  wheel  is  in  converging  beam  space.  When  this  is  combined  with  a  vari¬ 
able  rate  in  phase  wheel  rotation,  the  turbulence  can  be  fine  timed  to  emulate  a  large 
parameter  space  of  turbulent  conditions.  Additionally,  the  phase  wheels  themselves 
can  be  exchanged  for  an  even  greater  degree  of  flexibility. 

6.2  Methodology 

The  experiment  was  designed  to  replicate  the  scenario  used  in  the  wave-optics  sim¬ 
ulation.  To  accomplish  this,  several  key  parameters  were  used  to  scale  the  experiment 
to  the  simulation.  However,  the  limitations  of  the  lab  setup  required  some  deviations 
from  the  simulation.  These  key  parameters,  and  the  method  used  to  calculate  them 
from  the  experimental  data,  are  described  below.  After  that,  the  procedures  used  to 
capture  the  data  are  described. 
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Table  2.  A  comparison  of  the  key  parameters  used  in  the  simulation  versus  those 
used  in  the  lab  experiment.  Those  marked  with  a  *  are  the  non-dimensional  scaling 
parameters. 


Parameter 

Simulation 

Experiment 

Transmitter  Gaussian  Beam  Waist,  W0 

5.0cm 

0.45mm 

Wavelength,  A 

1550nm 

1550nm 

Path  Length,  L 

100km 

7.78m 

Fresnel  Ratio*,  A0 

19.7 

19.0 

Rytov  Variance*,  of 

0.92 

2.6 

D/r0* 

0.32 

0.32 

6.2.1  Key  parameters. 

A  comparison  of  parameters  used  in  the  simulation  and  those  used  in  the  lab 
experiment  are  shown  in  Table  2.  The  turbulence  parameters  of  and  D/tq  are  actual 
measured  values  for  the  experimental  column.  For  the  experiment  and  simulation  to 
model  the  same  scenario,  it  is  required  that  the  Fresnel  ratio,  Rytov  variance,  and 
D/tq  are  the  same. 

Determining  the  values  for  the  experiment  were  straightforward  for  the  wavelength 
and  the  transmitter  beam  waist.  The  wavelength  of  the  Thorlabs  LPS-PM1550-FC 
is  1550nm,  and  the  output  of  the  Thorlabs  F230FC-1550  fiber  collimator  is  fixed  at 
a  diameter  of  0.9mm.  However,  due  to  the  compressed  design  of  the  ATS,  the  path 
length  is  not  just  a  simple  measure  of  the  distance  from  the  output  of  the  collimating 
lens  to  image  plane  of  the  imaging  system.  The  method  used  for  determining  the 
path  length  takes  advantage  of  the  Gaussian  beam  propagation  equations  and  solves 
for  the  propagation  distance.  Recalling  Eqs.  (23)  and  (26)  the  effective  propagation 
distance  Le  can  be  solved  for  using  the  measured  propagated  beam  waist  via 


Le 


kW0^W2  -  W02 
2 


(144) 


The  beam  waist  was  determined  by  fitting  a  Gaussian  profile  to  the  received  intensi- 
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ties.  Once  the  effective  propagation  distance  was  determined,  the  Fresnel  ratio  was 
calculated.  Adjustments  were  then  made  to  the  propagation  path  in  an  attempt  to 
get  the  Fresnel  ratio  as  close  to  the  wave-optics  simulation  as  possible. 

Determining  the  turbulence  parameters  required  measuring  the  resulting  beam 
waists  after  propagating  a  beam  through  the  turbulence.  Calculating  the  atmospheric 
coherence  length  r0  starts  with  the  relationship  between  the  long-term  propagated 
beam  waist  Wlt  and  the  short-term  beam  waist  Wst  given  as 


Kt  =  fVsr  +  ( rl ) , 


(145) 


where  (r^)  is  the  second  moment  of  the  instantaneous  center  of  the  beam  in  the 
receiver  plane  [1].  The  beam  wander  can  be  expressed  in  terms  of  ro  via 


\/(rf)  =  0.69L 
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(146) 


Using  the  effective  propagation  distance  in  place  of  L  and  solving  for  the  experimental 
atmospheric  coherence  length  r0  e  yields 


ro,< 


=  2IDn 


VWIt-W2st2W0 


-6/5 


0.69LP 


(147) 


Since  the  experimental  aperture  was  implemented  digitally  using  a  mask,  the  value 
for  experimental  D/tq  was  made  to  match  the  simulation  exactly  by  making  the 
experimental  aperture  size 

De  =  r0)e^-.  (148) 

The  other  turbulence  parameter  used  to  describe  the  experimental  configuration  is 
the  plane- wave  Rytov  variance  a\  defined  in  Eq.  (55).  Using  Eq.  (91)  and  assuming 
the  transmitter  is  fixed  and  assumed  to  have  no  mechanically  induced  jitter,  the 
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effective  long-term  beam  waist  can  be  calculated.  Then  solving  for  erf  yields 


(149) 


where  W  is  the  vacuum- propagated  beam  waist  given  by  Eq.  (26). 

6.2.2  Experimental  procedure. 

The  execution  of  the  experiment  consisted  of  the  initial  configuration,  the  data 
collection,  and  analyzing  the  data.  The  initial  configuration  consisted  of  aligning  the 
beam  and  determining  the  effective  propagation  distance.  The  Thorlabs  F230FC- 
1550  fiber  collimator  was  placed  in  a  Thorlabs  K6X  six-axis  kinematic  optic  mount. 
The  K6X  was  mounted  on  an  elevation  stage,  which  was  then  mounted  on  top  of 
two  translation  stages  allowing  for  a  three-dimensional  translation  adjustment.  The 
K6X  allowed  for  small  translations,  but  not  large  enough  for  the  experiment.  Once 
the  transmitter  was  mounted  and  aligned,  the  measurement  of  effective  propagation 
distance  was  made  by  capturing  vacuum  propagation  images  and  using  Eq.  (144). 
Adjustments  were  then  made  to  the  location  of  the  imaging  system  until  the  effective 
propagation  distance  was  8.1m,  the  distance  needed  to  produce  the  same  Fresnel 
ratio  as  the  wave-optics  simulation.  However,  the  technique  used  to  measure  the 
propagation  distance  appeared  to  be  very  sensitive  and  only  repeatable  to  about  95% 
accuracy.  The  value  shown  in  Table  2  was  the  measured  value  obtained  during  the 
data  collection. 

Once  the  initial  configuration  was  complete,  the  data  taking  process  began.  To 
take  advantage  of  the  same  routines  used  to  analyze  the  wave-optics  simulation  data, 
the  same  seventeen  beam  locations  were  used  in  the  experiment.  The  beam  separation 
used  was  2 pc,  the  same  as  the  simulation.  This  resulted  in  an  experimental  distance 
of  1.4mm  which,  when  considering  the  factor  of  four  demagnification,  translated  into 
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Figure  27.  The  mount  used  for  the  transmitted  beam  consisting  of  an  elevation  stage 
on  top  of  two  translation  stages.  The  configuration  allowed  for  translation  adjustments 
in  three  dimensions  and  tip-tilt  adjustments  in  two  dimensions. 


23.6  pixels.  A  table  was  then  created,  reproduced  in  Table  3  of  the  positions  needed 
to  capture  all  seventeen  locations.  Although  the  mount  was  very  stable  once  moved, 
moving  the  beam  in  a  repeatable  manner  proved  to  be  very  difficult.  Eventually  a 
technique  where  every  movement  of  the  beam  was  made  such  that  the  last  turn  of 
the  knobs  were  always  in  the  same  direction.  Although  the  angle  of  the  mount  would 
change  when  the  knobs  were  turned,  this  technique  helped  to  keep  the  final  angle 
the  same  after  each  movement.  The  beam  position  was  verified  by  moving  the  phase 
wheels  of  the  ATS  in  their  vacuum-propagation  position  and  observing  the  centroid 
shown  on  the  camera  interface.  Once  the  beam  position  was  verified,  the  phase  wheels 
were  put  back  into  their  turbulence  position.  To  collect  the  data,  one  phase  wheel 
was  moved  at  0.5RPM  and  the  other  phase  wheel  was  moved  at  5.0RPM.  Camera 
images  where  saved  over  a  period  of  125s,  ensuring  that  a  full  rotation  of  the  slower 
phase  wheel,  and  twenty  rotations  of  the  faster  one,  were  captured.  In  addition 
to  the  seventeen  beam  positions,  a  vacuum-propagated  beam  was  captured  at  the 
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Table  3.  Table  showing  the  relative  and  absolute  locations  of  the  beam  center  at  the 
receiver.  The  table  assumes  the  detector  is  centered  on  pixel  158, 125. 


Position 

Relative 

Absolute 

X 

Y 

X 

Y 

1 

0.0 

0.0 

158.0 

125.0 

2 

-11.8 

0.0 

146.2 

125.0 

3 

11.8 

0.0 

169.8 

125.0 

4 

00.0 

13.6 

158.0 

138.6 

5 

-11.8 

-6.8 

146.2 

118.2 

6 

11.8 

-6.8 

169.8 

118.2 

7 

0.0 

20.4 

158.0 

145.4 

8 

0.0 

-20.4 

158.0 

104.6 

9 

-23.6 

13.6 

134.4 

138.6 

10 

23.6 

13.6 

181.6 

138.6 

11 

0.0 

-27.2 

158.0 

97.8 

12 

-11.8 

20.4 

146.2 

145.4 

13 

11.8 

20.4 

169.8 

145.4 

14 

-23.6 

0.0 

134.4 

125.0 

15 

23.6 

0.0 

181.6 

125.0 

16 

-11.8 

-20.4 

146.2 

104.6 

17 

11.8 

-20.4 

169.8 

104.6 

center  beam  position  as  well  as  a  background  image  with  the  laser  off.  The  averaged 
background  image  was  subtracted  from  each  of  the  beam  positions,  including  the 
vacuum-propagation,  prior  to  analyzing  the  data. 

The  data  were  analyzed  in  the  same  manner  as  the  wave-optics  simulation  data. 
After  determining  the  proper  receiver  diameter  to  match  the  simulation’s  D/tq,  a 
mask  was  applied  to  the  images.  The  irradiance  for  each  frame  was  calculated  by 
spatially  integrating  the  pixel  intensities  across  the  masked  images.  This  process 
resulted  in  an  irradiance  matrix,  which  was  then  multiplied  by  the  evaluation  matrix 
in  Eq.  (135)  to  produce  the  irradiances  for  one  to  seven  beams. 
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Received  Irradiance  [dB] 


Figure  28.  CDF  of  integrated  irradiances  for  one  to  seven  beams.  Power  is  shown  in 
dB’s  below  the  mean  irradiance  for  a  single  beam. 

6.3  Results 

The  results  of  the  experiment  for  all  seven  configurations  are  shown  in  Fig.  28.  The 
results  are  shown  as  a  CDF  normalized  to  the  mean  irradiance  for  a  single  beam.  By 
normalizing  in  this  way,  the  shifts  in  the  CDF  due  to  changes  in  the  mean  irradiance 
from  the  different  beam  configurations  is  preserved. 

Figure  29  shows  a  comparison  of  the  multi-beam  model  with  the  experimental 
data.  An  effective  C2  value  was  originally  determined  by  calculating  it  from  the  atmo¬ 
spheric  coherence  length  and  the  Rytov  variance  via  Eqs.  (43)  and  (55),  respectively. 
The  effective  constant  C2  calculated  from  the  atmospheric  length  was  6.8  x  10-11m-2/3 
while  the  value  calculated  from  the  Rytov  variance  was  9.8  x  10_1°m-2/3.  However, 
the  value  determined  to  fit  the  single  beam  CDF  best  was  C2  =  9.0  x  10~nm~2/3. 
Both  methods  of  estimating  the  turbulence  are  very  sensitive  to  the  estimate  of  the 
effective  beam  waist  and  the  initial  collimated  beam  waist.  However,  by  using  the 
single  beam  CDF  to  estimate  the  value  of  C2  a  good  comparison  can  be  made  for  the 
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Figure  29.  CDF  of  integrated  irradiances  for  one,  two,  four  and  seven  beams  compared 
to  model.  Power  is  shown  in  dB’s  below  the  mean  irradiance  for  a  single  beam. 

multiple-beam  cases. 

6.4  Discussion 

The  model  produces  a  reasonable  prediction  of  the  experimental  results.  Differ¬ 
ences  can  partly  be  attributed  to  some  of  the  unavoidable  inconsistencies  in  the  lab 
setup.  Particularly,  the  method  for  positioning  the  beams  may  induce  some  point¬ 
ing  error  resulting  in  the  beams  not  traveling  in  parallel.  Additionally,  due  to  the 
nature  of  the  stepper  motors  used  on  the  phase  wheels,  the  slow  rotation  of  the  first 
wheel  caused  a  small  amount  of  jitter  in  wheel’s  rotational  rate.  Finally,  the  camera 
frame  rate  was  inconsistent  due  to  the  lack  of  an  external  trigger.  This  variation  in 
the  frame  rate  would  cause  the  post-processed  integrated  irradiances  to  not  properly 
average  out  the  turbulence. 

The  model  also  contains  some  assumptions  that  are  not  present  in  the  experi¬ 
ment  or  the  simulation.  The  model  essentially  assumes  that  each  transmitter  is  the 
same  radial  distance  from  the  receiver,  which  for  more  than  three  transmitters  is  not 
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true.  It  does  this  by  weighting  their  contribution  to  the  new  shaping  parameters  in 
Eqs.  (127)  and  (128)  equally.  The  wave-optics  simulation  violates  this  assumption  in 
the  same  way,  and  experienced  a  similar  disagreement  with  the  model  as  the  number 
of  transmitters  increased. 

In  summary,  the  performance  of  the  new  model  has  been  assessed  by  comparing  it 
with  results  from  a  wave-optics  simulation  and  a  lab  experiment.  While  some  of  the 
assumptions  used  by  the  model  are  not  valid  in  real  world  scenarios,  the  model  still 
provides  a  useful  tool  for  predicting  the  performance  of  a  multi-beam  FSOC  system 
in  turbulence.  Additional  testing  of  the  model  would  be  beneficial,  particularly  in 
a  larger  held  experiment  which  would  provide  better  scaling.  A  held  experiment 
would  also  be  able  to  use  multiple  transmitters  at  the  same  time,  which  would  hx 
the  alignment  issues  associated  with  moving  a  single  beam  to  the  different  beam 
locations. 
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VII.  Conclusion 


This  chapter  discusses  the  final  research  conclusions.  It  includes  some  of  the 
challenges  that  were  overcome,  the  key  research  results  and  contributions,  and  some 
recommendations  for  future  work. 

7.1  Challenges  overcome 

The  first  challenge  involved  defining  the  problem  to  be  researched  and  scoping 
the  problem  in  such  a  way  that  the  research  would  be  completed  in  the  time  avail¬ 
able.  Building  on  previous  free  space  optical  communication  (FSOC)  work  on  from 
recent  years,  some  fruitful  research  areas  were  found  where  overall  performance  im¬ 
provements  could  be  achieved.  Previous  multi-beam  work  was  limited  to  using  two 
transmitters,  with  the  use  of  a  larger  number  of  transmitters  and  their  configuration 
remaining  uninvestigated. 

The  next  challenge  resulted  from  the  lack  of  a  verified  probability  distribution 
function  (PDF)  model  for  the  irradiance  resulting  from  multiple  Gaussian  beams. 
Fortunately,  the  development  for  the  sum  of  independent  Gamma-Gamma  random 
variables  was  accomplished  by  Chatzidiamantis  [8]  shortly  after  this  research  be¬ 
gan.  By  applying  the  foundations  of  single-beam  irradiance  PDFs  from  Andrews  and 
Phillips  [1,2]  and  assuming  the  beams  were  independent,  the  combination  provided 
a  bounding  solution.  The  resulting  model  was  verified  with  wave-optics  simulations 
in  the  course  of  this  research.  This  model  provided  the  basic  framework  for  partially 
correlated  beams.  The  method  of  determining  the  correlation  coefficient  as  a  func¬ 
tion  of  beam  separation  provided  the  necessary  foundation  to  determine  the  partially 
correlated  PDF.  Finally,  weighting  the  contribution  of  the  averaging  effect  between 
the  beams  by  a  function  of  the  correlation  coefficient  ultimately  resulted  in  a  new 


84 


model,  developed  under  this  research,  for  multiple,  spatially  separated  beams. 

Another  challenge  involved  settling  on  a  method  for  propagating  multiple  beams 
in  a  wave-optics  simulation,  and  the  beam  configuration.  The  orientation  for  two  or 
three  beams  is  straightforward.  However,  the  use  of  four  or  more  beams  required 
a  method  to  ensure  the  necessary  beam  spacing  was  maintained  while  minimizing 
the  physical  footprint  required  to  house  all  of  the  beams.  By  using  a  hexagonal  grid 
structure,  the  effective  footprint  of  any  number  of  beams  is  minimized.  The  beam 
positions  are  also  such  that  the  beam  configurations  are  centered  on  the  receiver  and 
preferably  symmetrical.  These  constraints  led  to  the  configurations  shown  in  Fig.  13. 
To  reduce  the  computation  time  required,  the  union  of  beam  positions  among  the 
seven  different  configurations  was  determined  resulting  in  the  beam  positions  shown 
in  Fig.  18.  This  ultimately  reduced  the  number  of  individual  beam  propagations  from 
28  to  17  reducing  computation  time  by  40%. 

The  last  big  challenge  was  determining  methodology  for  propagating  multiple 
beams  in  the  lab  experiment.  The  close  beam  spacing  in  the  scaled  experiment  pre¬ 
vented  them  from  simply  being  placed  next  to  each  other  on  a  mount.  Initially,  a 
transmitter  configuration  was  developed  where  up  to  four  beams  could  be  propagated 
simultaneously.  However,  maintaining  beam  alignment  and  spacing  proved  to  be  ex¬ 
tremely  difficult.  This  resulted  in  using  a  single  beam  which  was  iteratively  moved  to 
the  different  beam  positions,  similar  to  what  was  done  in  the  wave-optics  simulations. 
The  initial  mount  did  not  provide  the  accuracy  needed  and  was  too  limited  in  trans¬ 
lation  distance.  Placing  the  mount  on  translation  and  elevation  stages  solved  this 
problem,  but  introduced  additional  pointing  errors  that  required  some  attention  to 
ensure  the  error  could  be  minimized  and  the  beam  positions  were  repeatable.  Better 
fiber  launch  mounts  may  have  eliminated  some  of  the  alignment  problems,  but  such 
systems  were  well  outside  the  available  budget. 
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7.2  Key  research  results 


This  section  summarizes  the  key  results  and  contributions  of  this  research. 

7.2.1  Analytical  PDF  for  multiple  beams. 

This  work  provides  the  first  development  of  an  analytical  PDF  model  for  a  multi¬ 
beam  FSOC  system.  The  initial  model  assumes  the  beams  are  completely  indepen¬ 
dent,  providing  a  best-case  scenario  that  bounds  the  performance  of  a  multi-beam 
system.  Building  on  previous  work  [31-33,51],  the  necessary  beam  separation  was 
then  determined.  This  enabled  development  of  a  PDF  model  for  multiple,  spatially 
separated  beams.  Ultimately,  this  provides  an  engineering  tool  for  performing  design 
trade-offs  between  the  number  of  beams  used  and  other  FSOC  performance  parame¬ 
ters. 

7.2.2  Improved  FSOC  performance  with  PPM. 

This  work  demonstrated  FSOC  BER  improvement  using  pulse  position  modula¬ 
tion  (PPM)  for  FSOC  when  experiencing  scintillation.  The  obvious  performance  im¬ 
provement  is  seen  in  Figs.  20  and  21  which  demonstrate  the  utility  of  using  PPM.  Pre¬ 
vious  techniques  using  adaptive  thresholding  [33]  with  on/off  keying  (OOK)  showed 
some  improved  performance.  However,  PPM  inherently  provides  results  that  are  con¬ 
sistent  with  what  would  be  achieved  using  OOK  with  instantaneous,  perfect  adaptive 
thresholding.  The  cost  of  this  improvement  is  a  drop  in  bandwidth  efficiency  by 
a  factor  of  two.  Additionally,  by  combining  multiple  beams  with  PPM,  significant 
BER  performance  gains  can  be  achieved.  This  integration  addresses  the  two  aspects 
of  scintillation  in  FSOC,  including  reducing  the  probability  of  fades  and  mitigating 
the  impact  of  fades  on  overall  BER. 


7.2.3  Verification  with  experimental  results. 


This  work  provides  experimental  validation  of  a  multi-beam  model  for  FSOC  ap¬ 
plications.  Previous  work  predominantly  relied  on  on  wave-optics  simulations  to  ver¬ 
ify  multiple  beam  propagation  effects  [31-33,51].  When  used,  experimentation  with 
multiple  beams  was  generally  done  for  other  applications  [23]  and  not  FSOC.  Ex¬ 
perimental  validation  in  the  context  of  FSOC  systems  further  solidifies  the  improved 
irradiance  statistics,  namely  the  reduced  scintillation  as  the  number  of  transmitters 
is  increased. 

7.3  Recommendations  for  future  work 

During  the  course  of  this  research  several  areas  for  future  research  were  identified. 
These  areas  are  summarized  below. 

7.3.1  Transmitter  tracking  and  pointing. 

Both  the  simulation  and  lab-based  experimentation  in  this  work  was  based  on  fixed 
transmitters  and  receivers,  with  transmitted  beams  being  parallel  to  each  other.  The 
lab  experiments  highlighted  the  practical  difficulty  with  aligning  multiple  beams  such 
that  they  travel  along  parallel  paths.  In  practice,  a  tracking  system  at  the  transmitter 
end  would  need  to  track  the  receiver  to  ensure  proper  beam  pointing.  This  could  aid  in 
keeping  the  beams  parallel,  possibly  by  calculating  the  required  offset  needed  for  each 
beam.  Another  method  would  be  to  simply  allow  the  beams  to  point  to  the  receiver, 
resulting  in  converging  rather  than  parallel  beams.  This  would  require  a  tracker  for 
each  transmitted  beam.  Regardless  of  the  technique  chosen,  the  model  would  need  to 
be  updated  to  account  for  not  only  the  increased  mean  irradiance,  but  for  the  change 
in  scintillation  due  to  the  pointing  and  tracking  process  itself.  The  tracker  would 
effectively  track,  and  therefore  point,  bright  spots  in  the  scintillation  pattern  which 
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might  have  the  overall  impact  of  reducing  scintillation  at  the  receiver  [14]. 

7.3.2  Receiver  optimization. 

The  receiver  used  for  wave-optics  simulation  and  lab-based  experimentation  used 
a  direct-detection  technique.  While  this  allowed  the  research  to  concentrate  on  the 
irradiance  statistics  of  multiple  beams,  it  did  not  account  for  phase  perturbations  at 
the  receiver.  If  the  receiver  were  to  couple  the  light  onto  a  fiber,  or  focus  the  light 
to  some  other  free-space  detector,  significant  signal  degradation  would  be  introduced 
based  on  wavefront  distortions. 

In  the  case  of  multi-beam  FSOC,  the  receiver  design  would  also  have  to  account 
for  the  received  image  of  the  multiple  beams.  If  the  receiver  were  able  to  resolve 
the  different  beams  in  the  far-held,  only  light  from  the  center  beam  may  be  coupled 
into  the  detector.  Additionally,  if  the  transmitters  were  designed  to  point  towards 
the  center  of  the  receiver,  each  beam  would  reach  the  detector  with  different  phase 
fronts.  This  will  generally  impact  performance  of  a  fiber  coupler. 
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